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Sists of a combined Closed Brayton cycle and bottoming Rankine 
cycle. Coal combustion capability is achieved by using a 
Pressurized Fluidized Bed (PFB) heater. A thermodynamic 
analysis and volume estimate was conducted for various Brayton 
cycle compressor pressure ratios (7, ). 


The results show that a maximum thermodynamic 
efficiency ( ?+ ) is achieved in the range of Tle = 5 to 6. 
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The design point was ultimately chosen at T. = 4.5 with a De = 
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technology is expected te come on line within the next few years. 
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iv NRO UCT LON 


ie Furpose. Increasing prices and decreasing supplies of 
fuel oil nas created a need for highly efficient ship pro- 
pulsion systems capable of utilizing fuels other than oil. 
The purpose of this thesis is to perform a preliminary design 
of a combined cycle propulSion plant equipped with a multi- 


fuels combustor capable of burning coal or oil. 


”~ 


2. system Description. The system investigated 1S a closed 
regenerative Brayton Cycle combined with a bottoming Rankine 
cycle. Working fluids for the cycles are air and water 
respectively. The heater for the Brayton cycle 1S a pressur- 

ized fluidized bed (PFB). The heater cycle is essentially an 

open regenerative Brayton cycle. This thesis was conducted 
considering the heater only in its coal burning mode of opera- 
tion. Work done by the Brayton and Rankine cycles will drive 

a common shaft through a double reduction type gear. An enthalpy- 


enthropy diagram is shown on Figure 1. 


oye Scope. 
ae QObdjectives. The objectives of this thesis are to: 
(1) Perform a thermodynamic analysis of the system 


for various closed Brayton cycle pressure ratios utilizing 





pisure 1. 


Heater Cycle 


Closed Brayton Cycle 








realistic values for component efficiencies and losses. 
(2) Estimate major component volumes and thus obtain 
a total volume for set of parameters listed in (1) above. 
(3) Select a set of parameters for design based 
upon the above results. 
(4) Discuss the technical feasibility and implications 


of the design aS a means for naval ship propulsion. 


b. Design Criteria. 


(1) The system parameters selected for design were 
those which gave the lowest total volume within an acceptable 
range of thermal efficiency. 

(2) An acceptable range of thermal efficiency is 


defined as 11% of the maximum efficiency. 


Se Limitations to Scope. 

(1) No component or system testing was carried out. 

(2) The system was designed with the aid of a 
computer utilizing performance parameters of currently avail- 
able equipment. 

(3) The concept of pressurized fluidized bed (PFB) 
combustion has not been proven in a shipboard environment. 
The thesis was conducted based on the assumption that a PFB 


System would be put to sea in the near future. 





II. THERMODYNAMIC ANALYSIS 


me rne Brayton Cycle 
a. Working cycle. This is a closed Brayton cycle 


utilizing air as the working fluid. Air was chosen since 
this is proposed for surface ship application and would thus 
eliminate the need for carrying a reserve supply of gas if 
something like Helium or CO. were used. 

Advantages of a closed cycle over an open cycle are 
the following: 

(1) Compressor suction pressure is no longer limited 
to ambient. By increasing the working pressure of the cycle 
a Significant reduction in engine size can be achieved for a 
specified output. 

(2) Higher working vcressures also increases the heat 
transfer capacity of the working fluid thus resulting in 
Smaller, more effective heat exchangers. 

(3) The use of a high pressure fluid gives the added 
advantage of regulating the cycle pressure level by varying 
the inventory of gas in the loop. If this is done in a manner 
which kept the maximum temperature and speed constant, the 
turbomachinery vector diagrams would remain essentially un- 
changed. The result is little or no change in compressor or 
turbine efficiencies with varying load, thus maintaining high 
part load efficiency. 

A Temperature - entropy (T-s) diagram of the working 


cycle is shown in Figure 2. 








Figure 2. 


point 1: compressor suction, cooler discharge 
2: compressor discharge, regen inlet (cold) 
3: regen outlet (cold), heater inlet 
4; heater outlet, turbine inlet 
Be turbine outlet, regen inlet (hot) 
6: regen outlet (hot), WHB inlet 
fe NHB outlet, cooler inlet 
Qn? heat input from heater cycle 
Qn? heat released to Rankine cycle 
QQ: heat released to sea water cooler 
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The analysis was conducted by assuming values for 
Ty» Ty Ne (p)( compressor polytropic efficiency), N+ (P) (turbine 
polytropic efficiency), En (regenerator effectiveness), and 


).. (total pressure drop). 


The specific heat was assumed to vary with temperature 
according to the following relation: 
Cp(T) = .2475 - (3.759x1077) T + (5.106x107°)P*-(1.231x107!4 )r? 


for T in OR 


Cp in BTU/lbm-°R 


The following are defined: 


oi - compressor pressure ratio (P/P, ) 
fim turbine pressure ratio (Pe/Py,) 
(Af) - % total pressure drop around closed Brayton loop 
Wp - power (net) from closed Brayton cycle 
Mi, - mass flow of working fluid in Brayton cycle 
A - temp in °R at point i 
ns ~ enthalpy at point i 
Ig - gas constant 
P. - pressure abs. at point i 
a - Specific volume 


ER - regenerator effectiveness 


The basic procedure was to determine To» Bas Ts» 6? 


Q;,' and W, for various pressure ratios (TT). 


ii 





Determination of T, cae icon) : 


Small stage polytropic efficiency is 


wd 


~ Clie 
De (P) ~ ae Cp dT 


Since for an isentropic process nae = O=-dh;.- vdp 


Suostituting VU = RT/p gives 


maus, 








ie 7e 


Te ; 
So es eee ea a 
fs Se dk a BZ ie ~ 2p) fn te 


y 


By substituting the equation for or C(T) one gets 





2675 w(z2\; 3.759x10-9(7,-T,) + (2.553x107°)(1,°-7,°) - 
T, 
1 
(4.103x107°*) (757 - 7,7) 
ves An The 
/<(P) 


Ts Can then be determined iteratively. 


Determination of bgt 


In a manner Similar to the compressor above, 


Rs 


5 
§ 2. AT 7 R ap An (Be 


= 
= 12 








After integrating, this becomes 


2 va 
-1),“)- 


~5 HO 
2475 (es) (3.759x10°°)(T,.-T, )+(2.553x10°-~)(T, 
4 


(4103x107 *) (0,7 - Ty?) 


= RQ ln(Mty) 


T. can then be determined iteratively. 


MeLtermination of Us and Te: 


i and T¢ can be found from an assumed regenerator 


3 


effectiveness (Ep) using the following relation: 


Bae Oe aes : Ale or 
ts - T 5 T. - T. 
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_ 7 
a4 
"2399 akat | | 
7 ~~ 


"0 DeRaGh a fs 
. 


b. Heater Cycle. This is an open Brayton cycle also using 
combustion air as the working fluid. The heater itself is 
a pressurized fluidized bed utilizing coal as fuel. A T-s 


diagram of the heater cycle is shown in Figure 3. 


Major components of the cycle consist of a heater, re- 
cuperator and turbine driven compressor. The turbine is 


powered by exhaust gases leaving the recuperator. 


The analysis was conducted by assuming values of Tah 
Toy? VelP), D4(P), and @ (equivalence ratio). Both the 
exhaust products and combustion air were assumed to be perfect 
gases with specific heats varying in accordance with the re- 


Mition given in i.(a). 


The following are defined: 





Ten - heater compressor pressure ratio P,,/P,,) 
TE - heater turbine pressure ratio (Pop /Pey) 
OP 

cn - % total pressure drop around heater loop 
me - mass flow of air in heater cycle 
Me - mass flow of fuel in heater cycle 

(F/A)> - stoichiometric fuel - air ratio of fuel used 
LHV - Lower Heating Value of fuel used 
Ty - absolute temp at point i, 
As - enthalpy at point oy 

? H - efficiency of the neater cycle 

Barc - recuperator effectiveness 
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ambient 





Rareate 3 « 


Meant Ih: compressor inlet 
ech: compressor outlet, recup inlet (cold) 
3n: recup outiet (cold), heater inlet 
Hh: max air temp in heater 
on: heater outlet, recup inlet (hot) 
6h: recup outlet (hot), turbine inlet 
fa: turbine exhaust 


3! heat released to working cycle 


an! heat in = Mp LHV 


15 





Assume the heater cycle to be a control volume (C.V.) 


as shown in Figure 4. 


mM p-LHV | : 
SSS, r——— 22 
| | 
Cvs 
| | 
| 
‘— & Oa | 
mM, Din (m, + ms) Don 
Figure 4. 


a, leaving the C.V. must be equal to Qn entering the 


working cycle thus, 
Qp = My (hy, - hs) alo + Mp) (hy - he) 
Performing a heat balance on the C.V. gives us 


ay 14 (m, + Ma) Don = m.(LHV) 33 moh. 


Therefore, 








= Gs = oe. = | ~[(1+(%)s 2) ho 5 hi] Al, OLAV] 


“) » 
Qin ne 
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From a power balance between the turbine and compressor 


1) An (no net work) 


ee be he = pa? = (my, + Oe) Com (Ten - Toy) 


put 
C Dea ) 
Tah) Salis (Tei) Pe /; be 
; ~1 
Web, Tee (Tn) ae Mh * (| [ I- Por [Treh) 
thus, 


Cp. Tin ten) 7 - || _ [| uy (A) D | Ce Th (fre) * i | 


=fi+ (A), BT Con, 4 EE) HOE Cet Ht K 


or, 


Che Tih rex) */ 1] +[ 1+ (FA). P| Co Tn = 


[1+lA), 9] CopTn [1-[F le 1M a GH 


ae 





therefore, 


AP See 
Dic ener) P| Coe Ta Fay — 


For minimum volume heat exchangers we want the largest 


A 
(“%) compatible with the above eqn. 


Gir rerentiating, 
“AP \ 
Con 1 F) Cpa -| [-3( Fm 
Al \ Oo Tek 


i{I+ (., o| Coy _ pipe ee Tip Co /' Ie ra io Sf (A), © 2] 


Coy Toph Coy ? (Ts,) = ay tie roan (Trek “ft oS || + { (ta ), ?|C Ces. Tra)$ 
i 1), Ge Ta Tay PPP? 


Setting (42) 
Th = 0 and solving for [ich gives 
Toh 


eee Hie lA), BJ See (Cu Tn - oy ST he [ep 


Che Tin (Sp. [Ne - Cog %) 


ue 





By substitution 


a. 7 dh fae 
far) =) - (Leta P Venta Tala] 9 
} [| uv TAY « D | Coy Tay (Tes) Ce/74 


ee Th 


AP 
Thus Net and (F-), can be found through known or assumed 


quantities. 


Determination of T 3h? on 


Si (Ts ~Téh oe ‘es GER - Tok 


rec 


Cmin (Teh - Teh) Comin (Toh - Teh) 
where Ch = Coy (mr. + mf) 
Ce a Cpe. uo. 
therefore, 


T3} = Teh + { + (ir) P| (Top GS) 
To, and Tg, can be found once Te, is determined. Assume 


that 


Ia 100°R 


5h 


By 





2. The Rankine Cycle. This is a bottoming steam cycle that 
recovers heat from the closed Brayten cycle that would ovther- 
Wise be discharged overboard through the cooler. The rela- 
tively hot gas leaving the regenerator is pre-cooled in the 
Waste Heat Boiler (WHB) before being sent to the cooler which 
brings the temperature down to compressor inlet conditions. 


A T-S diagram of the Rankine cycle is shown in Figure 5. 





eee! 5. 
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Major components of the cycle are the WHB, condenser, 
steam turbine, and feed pump. I% is assumed that the feed 


pump is driven with power from the turbine. 


The analysis was conducted by assuming values of Ty 5? 
aT, pls) (pump isentropic efficiency), b+ WS) (turbine isen- 


tropic efficiency), and y WHB (waste heat boiler efficiency). 


The following are defined: 


Ass - enthalpy at point is 
Bec = abSOluve eemp ay DOINt 1S 
Ee ~ pinch point on Brayton side 
SF - Saturation point on steam side 
pp - mass flow ratio (:,/m,) 
x - quality of steam exiting turbine 
—_ - absolute pressure at point is 
Wp - power (net) from the steam cycle 
Di - thermal efficiency of the combined cycle 
Ma - mass flow of steam in Rankine cycle 


Rankine cycle parameters were chosen for maximum work. 
The procedure was as follows: 
i Compute ae from the relation 
Tag = T¢ - AT 


O 
(2) Assume a value for Hoes a < oe < 662.0°F 


el 





-3) 


(4) 


(8) 


(9) 


B(i0) 


(11) 


(12) 


Compute Top ion thevmweration 


Top = Top -AT 


Mecvain F Rope Upc oon Steam tables. 


SAT’ 


Compute Ans from the relation 


yssNag- Ht Cis) (nag -Ryc Us) ) get hy ('S ) from steam tables 


fimeain xX (turbine exit), Aye from steam tables 


Compute Nos from the relation 


h 2stls) > ie 
Dp lis) 


Compute © from the relation 
Vas ts, Done (he ~ how) 


sa hss ~ Asp 


ees ist 


— 


Compute Wp from the relation 


We = mg B[ (hss Figs) ~(ha¢ -h)| 


Iterate with T (step 2) until Wp is maximum. 
ee 


Compute T., from the relation 


_ a 
Cp (T,)+T; = Cp (Tre) - Tepe - 2 (he “hie )/ June 


ealculate Nes fesem tChemee la tion 


‘cee ali 
a3 
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Iti. VOLUME ANALYSIS 


1. Turbomachinery. An analytical procedure for the minimum 
volume design of turbomachinery was developed by Lee (ref. 
10). The following equations and mathematical expressions 
summarize that portion of the procedure which are applicable 
to compressor and turbine designs in this thesis. Volume 
obtained by Lee's method neglects the volume of casing, inlet 
and outlet plenums, and bearing assemblys. 

a. Turbine. Free vortex blading and axial exit stages 


wereassumed throughout the analysis. 


The following are defined: 


S ae - factor of safety 

U, ~ max. allowable tip speed 

0. - flew coefficient at tip 

Vn - hub—cieohatio at Inlet 

Ve ~ MUb=cip ravao at exit 

Cy ~ clearance factor 

Vo - specific volume at inlet 

hy - max tip radius 

h - Poly ctropue «coecriicient 

K - isentropic exponent 

Pm - density of blade material 

x - Stress concentration factor at root 
a - an Dig wedi LOM tac tor 

o ~ Gouvection ftacuor for taper 
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Xy - Cicreenonmteaece 1St principle axis 
chord 

Ar - normalized moment of inertia 

Vol, - volume of turbine 


The volume was computed by assuming values for SF, U,, D.. 


Veo» Cry Ks Pm % + Ay Xp» Keg and 0. 


Cis Ny and Yn are computed from the following relations: 


Faget (me Vox /f Pe Ue 1 (1- Ye") 1) 


aaa 
K ~(K~-1) 940) 


! en 
ye = [i-m* (-m2)]" 


Compute Cy Coy Cas Cy, from the following: 
cy = aya cca) Ma (1- Ya”) Uy" ro Xi / Ky 
c,= 2(A-BR) /(SF+«) 


Ces ZI (1 YX) /(sF +) 


Ca = - OPm Ue /2 


where A and B are material dependent constants. The following 


Can then be calculated: 


=¢,[(i-Y)/G-y)] o; 
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pac. tie) + Cali-y’) 


i) 


ee as | - (o + y?) hs y) + 2M Nf +)~(-D C- Woe “y?) 


a 


DD = 2V(I-Y*) “(n-) Cg - 2 C4 


The axial chord for a stage is given by DL. ROD da 
large number of stages, summation of Db, for the whole turbine 


can be approximated by an integral: 
yg 


Vis 
2 a e 


Ves 


where 


lon) = (%)” [(on-> -DD-N) /a*] ay 


This integral can be evaluated numerically thus, 


Vole = a ra a l, 


Dd. Compressor. An approach Similar to that of the 
turbine was used for the axial compresSor. Axial exit stages 
wereassumed in the analysis. 

The volume was computed by assuming values for SF, Uys 
Dar Yes Cr. Keo Xs ea | Xr Ks, » 9, and Vy (blade material 
yield stress). 


£5 





Tye Ne and Ve are computed from the following relations 


a (tr Ve /L P+ Ue TF (1-Ye*)| ‘is 


ae 
rae -(K-1) /%e(p) 


Ve | . 
Ya =Ll- TM (i - ¥) 
Compute Cys Cos and C. from the following: 


es Ss 
C=. 2 (Var DUR re De® XK /Ka 


C2 * -Ofm Uy /2 


ee = 2 Vy [x + SF) 


Then calculate the following: 


N= CLOW) (Oe 7 DMO? + 722597) 
m = Cs rC,(-9) 


Dy =C,{- (t+ AY) [as re2sy ry) 
— (1-9) (i eV)? (De? +499") Oe? +. 72259) 
+ 98 Vi-V) Ge /( G2 +. 7228 V) 
- 4.445 (-Y) ry! (oe? + FTV) (BE + 7225)" | 
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As with the turbine define 
y 
Ce 
I. = dX (b},) 
7 


where 


din = EU)” [(pn-d -po-n) b*] dy 


This integral can be evaluated numerically thus 


Vol, = ZW ay ile 


2. Heat Exchangers. The volumes of the Regenerator, Re- 
Cuperator, Cooler, and Heater were computed using the metnod 
of minimum volume design developed by Lee (ref. 10). In 
essence Lee found that the optimal diameter ratio (D,/D) for 
these heat exchangers can be established as a simple relation 
hoethe pressure ratio (TT) of the cycle. The results of the 
optimally distributed pressure drops (a) as reported by 

Lee were also studied. Parametric equations relating pressure 
drop to pressure ratio were developed and used in this analysis. 
No consideration has been made on computing the volume of the 
shell, headers, or other appendages. Lee's relations were 
developed from basic heat transfer and pressure drop relation- 


ships. 
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a. Regenerator. The regenerator is assumed to be a 
counter flow shell and tube neat exchanger. High pressure 
is inside the tubes as is the case for current design prac- 


fices. 


The following are defined: 


(25) _ =  % pressure drop across the regenerator 
7 Keg (both steams) 
D 
ES) - diameter ratio for regenerator 
Req 
Ak - enthalphy difference across the interior 
i Reg stream of Reg. 
D - tube inside diameter 
AT - log mean temp. difference of not and cold 
streams 
sa - Pranclt number 
4 - density 
KM - dynamic viscosity 
k - chenmal conductivity 


af {Dh 
Values ror Reg and (> ~ Reg Can be estimated from 


the following relations: 
neg SS be 0284 TT (Boia m = +4, Ep = .88) 


D; = 
) neg 7 +557 + 09337, 
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Subscript 1 refers to quantities inside the tubes 
whereas subscript 2 refers to quantities outside the tubes. 


A can be computed as follows: 


: eet ac ls os 2f Da\* 
A rec =(™%,) Ve Yn Dr 6 Ani . (1+ (42), [ ane ‘ca) | 


Reg 


where 
me (Ao = | 


-i.Fi i . 
=e A, FR -; 


-,4/ | i? —a -( 4 


ea De AT 


The volume can then be computed as follows: 


: AP\ 0.44 
Volpog ~ Reg we 


b. Recuperator. The recuperator is similar in design 


to the regenerator only it refers to the heater cycle. 


The following are defined: 


a esc - % pressure drop across the recuperator 
(vboth streams) 


(2) nec - dvameverwracio for recuperator 
‘AP \ - % pressure drop across the heater side 
Pp Shir, (outside tubes) 
(BF) ore - % pressure drop across the PFB 
An -  Enthalpy difference across the interior 
‘Rec 


Stream of recup. 
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AP\ Di 
Values tor (> rec and(2)e 4, can be estimated from the 


following relations 
a 2 (45 | f ( aP\ AP \ 
(4 Ree \ P STh ries -(- “p> PFE 


Dh, = lana 
Po) Rec 65 


Subscripts 1 and 2 refer to quantities inside and outside 


the tubes respectively. AR a6 1s computed as follows: 


D PL GEN -3 4.41 
Nee <0) TS te Ahn” [ie (S Deel [+B] 


where Y and Y are as previously defined. 


Then compute the volume 


. _ Ae Onl 
Ecc et Rec 


eae COoler. The cooler is assumed to be a shell and tube 


heat exchanger with the gas inside the tubes and salt water 


Sutside., 


The following are defined: 


Ar 
(> Pa. - % pressure drop across the cooler (air side) 
oir = diameter ratio for cooler 
Ah, eilg = enthalpy difference across the interior 
stream of cooler 
“WHB - sp. vol. of working fluid at WHB exit 
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AP 
Fins - % pressure drop across the WHB (air side) 


(Se) - % pressure drop across the heater 
ae (inside tubes) 
Mp (P) - polytropic cooling water pump efficiency 
Vig - Sp. vol. of cooling water 


(Perr ana (2 
Values for\?r/Clir and\_p/Cir are estimated from the 


following relations: 
(AP) =m sod - (42) aS), 
es WHE EO Reg Ete) 


Din 
Ae, Sore Mina. 7 © 
Subscripts 1 and 2 refer to quantities inside and outside 


the tubes respectively. Acir is computed as follows: 


Ne CA) Te Be tng Ate i (Be re (Ba) "(1+ (2 ‘) ‘< 
where V , °° are as previously defined and 
A(Ym) CW) Code) 
C= (Hin) I) 90) (“Nae )( Ze) 
OK MY tn) Ye) 


= 75 


oo 





The volume is thus 


: aN ola 
Voloir ~ Saw Se 


d. Heater. The heater is assumed to be a shell and 
tube heat exchanger with the working gas in the tubes and 


combustion gas out. 


The following are defined: 





DK ee 
nt - diameter ratio for heater 


Alby Tate ~ enthaphy difference across the interior 
; stream of heater 


Values of (45) a and (Dh) are 
PUSS OG ae nary D htr 


estimated from the following rélations: 


fi 


af) ipa) /AP\ = 
a 88 Me [for (SF) =. 1] 


3 
{AP AP " ) 
UP hte , ~ Ss htr 2 fp > htr 


where 


iS 


fp ("4a .) py” (WD 


(2h) = 049, +.2i4 


Subscripts 1 and 2 refer to quantities inside and outside 


pe 





the tubes respectively. Compute Antr as follows: 


F ° 
a 


aor, htr 


Lot 


Ait. =("4) ¥, Me Ahint- * [ : SS). E ea ayia | 


where Y , ¥ , % . are as previously defined. 


Then compute Voll sy 


is Die = 
Yitr = al ee, 


e. Waste Heat Boiler. The boiler is assumed to be a 
once through cross-counterflow type, as shown in Figure 6. 
High pressure is inside the tubes thus saturation pressure 
must be greater than the working presSure of the closed 


mmayton cycle. 
fas ou 


4 


Water in 


steam out 


gas in 


Figure 6. 
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Boiler volume is estimated using heat transfer relations 
and volumes from similar U.S. Navy COGAS (Combined Gas Turbine 


and Steam) designs (ref. 11). 


The following are defined: 


Cos = ave. specific heat on steam side 
Cy - ave. specific heat on air side 
NTU = Number Transfer Units 

#Ps - number cof passes 

Arm = Heat transfer area 

U - overall heat transfer coefficient 
Vin - total volume 


From a DTNSRDC (David Taylor Naval Ship Research and 
Development Center) study (ref.11) the following relations 


are obtained: 


CR = CMIN = ia Crs ob Ls Cp 
CMAX ms Cp m. Ces 


ER = (une OR) one ~1) 


NTU =C*R) An hi Yeala [I a0 2a" *) Kc0 -E07%)] 6 


_ NTU-CMIN 
Aur * U 





For similar designs assume Aum o< Vn 





ee OTs 
thus Vy = CMIN, ; NTU, 
Viz CMIN,+ NTUs + 
Vy (i CMIN, NTU, 
CMIN, 1 NTU. 
The vern d Can be evaluated from known data 
CMIN,, NTU+y 


on previous designs to give 


VHB ioe smo NEU 


But it snould be noted that the above equation assumes 


JIE . P 
constant (42)... and boller power. If as before Vol “™ Pann 


then from heat transfer 
thus 
Vo x m, Ah, 


Correcting the above equation for power gives 
Vig = ©0017 CMIN * NTU +m, * Oh, 


Ap \ 7! 
AS stated previously Vol ™ =) ; 
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Correcting again for pressure drop yields 


y 


’ AP 
wHR = 100045 CMIN * NTU rae mahe (S—) (£43) 





CMIN in BTU/s - °R 


m. in lbm/s 


2B Wass BTU/ibm 


ff. Condenser. The condenser iS asSumed to be a shell 
and tube crossflow heat exchanger with salt water inside the 


tubes and steam outsice. 


The following are defined. 


D - inside tube diameter 

N - number of tubes 

L - length of tubes 

g - Gravitational acceleration constant 

AP. - pressure drop on the salt water side 

19 ~ fore ulon. vac Lor 
te - density of cooling water 

- - velocity of cooling water 

ee cross sect. free flow area on water side 
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flow rate of cooling water 


heat flow in condenser 


w , (x) specific heat of cooling water (C,,) 
temp. difference across the cooling water side 


Steam quality entering condenser 


enthapy difference between sat. liquid and 
Sat. steam 


mass flow of steam entering condenser 


From basic heat transfer relationships 


where 


therefore, 


where 


Aum = Tf DNL 


H 


4Ac 
Teo 





Car. Dp 
ae re Vee 


ue Z2Ac g AR 
A iit 7 # PeVve~ 
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thus 


Dee oe x 
Ar ae Ve dee, 


; x 
For similar designs assume Vir Arm 


thus 
a Picts a g OP. Neg (rhs x); x): 
‘7 (RRL [eee ™ msn 
~ | for similar designs assuming same A P 
therefore 


b4 


= rs (ra, x). 
(ms ep 


The quantity in brackets can be evaluated from known 
Seta of previous designs (ref. 11) 


to give 
: 3 
leeng = 27 +3 mx ft 
m. in lom/s 
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veo bomeN OFA 25,000 Hr PLANT 


A proposed design of a 25,000 HP coal-fired combined 
cycle was completed. Coal comoustion is accomplished using 
a pressurized fluidized bed (PFE). Design parameters were 
eventually chosen after investigating the effect of changing 
pressure ratios, varying (AF) in the Brayton cycle, varying 
Ep (regenerator effectiveness), different compressor inlet temp. 


(T,) and steam condensing temp. (Tio), and variations in Brayton 


working pressure (WP) and saturation pressure CE acia rr, 


ae The Fluidized Bed. The pressurized fluidized bed 
operates at a pressure of 3 to 10 atmospheres. Pressure is 
maintained by an exhaust driven compressor. For proper com- 
pustion incoming air must be preheated by the recuperator to 
at least 700°F. The bed material is primarily limestone which 
serves to capture sulfer in the coal and thus minimize SO. 
emissions. A bed temperature of 1700°F is generally used to 
keep NO, emissions at an acceptable level. 

The fluid bed itself is about 4 ft. deep. Preheated 
combustion air enters the bed via a perforated bed plate under- 
neath the bed. Pressure drop across the bed is assumed to be 
10%, 

The PFB is used to heat the working fluid in the closed 


Brayton cycle. The working fluid leaves the regenerator and 
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enters the tubes in the PFB. The working fluid is first 
heated by convection from the gas leaving the bed and is then 
ducted to tubes submerged in the fluid bed for final heating 
to the turbine inlet temp (T)). (Figure 7) 


exhaust 
iN 


working fluid in(T. 





Pua bed 








{ 
Soa in. 


| == WOrKIng fluid our 
aie) 


mee i 
-_-_ —_- — a Ly 
Semoustion ————-» p- 
ee (Typ) crid plate 
ba Sure. 7s 
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b. AsSumptions. 


The computer program was run with the 


following input values: 


Ve (P) 
V+ (P) 
a Ee 
Dp Us) 


Dr) 

Sa 
PST 

? NAB 


m 
4 


i 


D(P)h 


0h 
g 


Risiie' 

91 

88 

80 

.00 

mic 8 

94 

540.0 OR 

2060.0 OR 

540.0 OR 

+35.0 OR 

1350 Si 

200.0 psi (want at least 50 psi greater 

than WP) 

25000 HP 

535.0 OR 

760 °R (limited by acid corroSion in the 
turbine) 

Po 5 

89 


.90 (11% excess air) 


In addition it was assumed that scoop injection could be 


used for both the cooler and condenser at speeds greater than 


10 knots, therefore power requirements for salt water circulating 
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pumps were excluded in the calculations. 


c. Results. The computer output is provided in Appendix 
B. Note that the saturation pressure for pressure ratios of 
2 and 3 were below the 200 psia specification. This is due 
to the fact that the low pressure stream leaving the regener- 
ator is at a temperature below the saturation temp. at 200pSia. 
Replot of V+ vs. 11 . and Vol, 4.) vs. iT, is shown in Figure 
Ne 


Note that the volume increases Sharply above ue = 4 
due largely to the rapid growth in the waste heat boiler. 
Maximum Vt OCCUns at around lige = 6. To investigate the 
effect of or on as and Vol., the program was run for 
2), roe pOO, .OC, «1O0,*and..l2 between Vis of 4 and 9. 


The results plotted on Figure 9 show that significant 
gains in !/, can be realized with very little increase in V, 
for the range of 4-6 pressure ratio. The effect of changing 
E. (regenerator effectiveness) in this range with a (4) = 


.05 was further investigated. 


The results listed in Table 1 show no increase in effi- 
ciency as En is reduced. Furthermore, volume increases Sharply 
aS a greater percentage of the work is being done by the Rankine 
cycle as Le decreases. 

Thus for greatest efficiency with minimum volume, cycle operation 


Should be at the highest Ee practicable. 
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Figure 9. 
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To investigate the effect of changes in qT, (compressor 
inlet temp.) runs were made on 4 different temperatures and 
the results shown in Figure 10. As expected efficiency de- 
creased with increasing rT, with very little change in the range 
of interest 4-4.5 Te: Volume on the other hand decreased 
Slightly for ie less than 3.5 but increased significantly 
above that. Thus for best efficiency and least volume Ty 


should be kept as close to cooling water temperature as possible. 


Finally, changes in working pressure of the Brayton cycle 
were investigated. Decreases in working pressure may require 
additional volume of Brayton cycle components yet at the same 
time allow a decrease in Rankine cycle steam vressure (recall 
50°R pressure differential between boiler streams). This will 
in turn increase cycle efficiency yet further increase volume. 


These results are born out in Figure il. 


d. System Selection. In accordance with the design 
criteria the highest thermal efficiency considered is 0.487 
(Figure i1). Thus it is desired to select the plant with min- 


imum volume at ke = een 


Figure 11 shows that only those alternatives with Brayton 
cycle working pressures of 70psi and 100psi meet tne efficiency 
requirement. Of these two the 100psi case has minimum volume 


at desired efficiency. 


However, the problem does not end here. If one looks at 
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the steam quality leaving the turbine at the like = 5, WP = 
moeenpSi point it is of the order of 0.84. This is a little 
lower than normally desired for design. In fact, the steam 
quality of all alternatives thus far considered within the 
range of interest is low. If a steam quality of order .87 is 
desired the current design pt. can change in one of 3 ways 
aS follows: 


1) Raise steam condensing temp (Ty) 


2) Lower minimum saturation pressure oe Aer, 


3) or a combination of the above two 


Through trial and error it was found that case 1 gave the 
greatest benefit with the least penalty. Computer runs were 
made in the 4-5 ii. range with the steam condensing temp 
raised to 560°R for both the 70 and 100psi alternatives. The 
results show a substantial increase in quality, slight decrease 


in efficiency and very little change in total volume (Table 2). 


(psi) mn 


Working Pressure ig 7 it x Vol, (ft) 
70 is 480 864 948.3 
V8 508) 482 Bewail 1140.6 
100 S20 479 ~862 956.5 
Table 2. 


The 70 psi designs have better steam quality and + than 
@he best 100 psi case. In addition, at 4.5 fe the volume is 
lower. Between the two 70 psi cases above a trade-off exists 
between 9 s and Vol,. Assume that volume comparable to the case 
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of i = 5S, = 100 psi, Ts = 540°R discussed above (Figure 
11) is desired. If this is the case then the new design point 
is at a Brayton cycle working pressure of 70 psi, lie = 4,5, 
Tis = 560°R. The results of this design are summarized in the 


following computer printout: 
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BRAYTON CYCLE: PRESSURE RATIO = 4.50 
WORKING PRESSURE = 70.0 

TOTAL PRESSURE DROP = .05 

MASS FLOW (LBM/S)=206.93 


STATE TEMP(DEG R) SUMMARY: THERMAL EFF=0.480 
i 540.0 HEATER EFF=0.933 (PHI=0.90) 
2 876.0 SFC (LBM/HP -HR )=0.440 
3 1403.9 VOL TOTAL(CU FT)= 948.3 
mn 2060.0 POWER (HP) =25000.0 
5 1475.9 WORK FRAC BY STM CYCLE=0.101 
6 948.0 
PP 836.2 
p 809.2 
COMPONENT DPP DHD  VOL(CU FT) 
. REGEN POLO. 9 50 26.5 EFF=0.88 
- TURBINE 5.5 EFF(P=0.91 
COOLER .009 .490 Welsieat 
COMPRESSOR 29.7 EFF(P)=0.88 
HEATER 0006 
BOILER .030 EFF =0.94 


HEATER CYCLE: PRESSURE RATIO=4.79 
TOTAL PRESSURE DROP =.43 

AIR FLOW(LBM/S)= 36.85 

FUEL FLOW(LBM/S)= 3.06 


Sale TEMP(DEG R) 


1H 535.0 
2H 855).2 
3H 1440.5 
4H 4809.8 
DH 1503.9 
6H 945.0 
7H Zoom 
COMPONENT DPP ODHD VOL(CU FT) 

RECUP 5 ol She oO, 1.9 EFF=0.91 
TURBINE 0.7 EFF(P)=0.89 
HEATER 019 «430 31.4 

COMPRESSOR 10.8 EFF(P)=0.85 


RANKINE CYCLE: SATURATION PRESSURE = 120.0 
STEAM FLOW(LBM/S)= 5.69 
FLOW RATIO =0.0275 


SALE TEMP(DEG R) PRESSURE ( PST) 


15 5OOn0 0.9 
2a oO T2020 
Sig SOle.Z 120.0 
35 913.0 1207.0 
4S Solel ee, 0.9 QUAL=0.864 
COMPONENT VOL(CU FT) 
BOILER 507.9 
TURBINE OE Pe is )=0.85 
CONDENSER gt 2 
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Via ae CATION 


The application of this system to Navy ships has many 
implications. One of which is the fact that neither Closed 
Brayton nor combined cycles have been previously used. However, 
there is a considerable effort at this time to get COGAS 
(Combined Gas Turbine and Steam) systems on board Navy ships. 
These systems are all limited to burning oil fuels. The Closed 
Brayton combined cycle when equipped with a multi-fuels heater 
offers the additional flexibility of burning the cheapest or 


most readily available fuel. 


Advantages of the Closed Brayton combined cycle are: 

1) The Brayton cycle operates at temperatures substan- 
tially below the material limit than in the open cycle gas 
turbines resulting in cheaper, more reliable, and longer life 
components. 

2) Low pressure steam systems are well proven, and have 
Significantly less water chemistry problems than higher pressure 
ones. 

3) AS load decreases and fuel flow is cut back, the con- 
trol valves (Figure 12) between the compressor outlet and cooler 
inlet can be opened as necessary (keeping max. temp and speed 
Semetant) So as to ensure that the vector diagrams of all sec- 
tions of the turbomachinery blading remain unchanged. Turbine 
and compressor efficiencies will remain essentially the same 


producing high ) at part load. This requires use of a CRP 
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Figure 12. 
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(controllable reversible pitch) propeller to maintain constant 
shaft speed at part loads. 

Disadvantages of the Closed Brayton combined cycle are: 

1) Higher overall system weight and volume due to the 
additional heat exchangers and piping required for the Closed 
Beayton cycle. 

2) A complex, integrated control system is required to 
coordinate control of the combined cycle, fuel flow, and CRP 
mo part loads. 

3) Fluidized bed combustion is not a fully developed 
engineering concept. Land based prototypes nave been built 
and operated witn limited but increasing levels of success. A 
shipboard prototype is perhaps 5 years in the future with pro- 
duction of a marinized version at least 10 years away. 

4) If coal is to be used as the primary fuel, it should 
be noted that it has a lower heating value per ton than oil. 
The result is a larger potentially more expensive (greater 
acquisition cost) ship for the same military mission. 

5) Reaction time with coal fired fluidized bed combustors 
to load changes are generally slower than that for oil heaters 
necessitating in a separate fuel oil system for rapid increases 


iy load. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The Closed Brayton combined cycle offers an attractive 
alternative system for many shipboard applications. As with 
the first Arab oil embargo, subsequent cut-offs of Arab oil 
could result in decreased steaming days and a lower overall 
readiness of the U.S. Navy's defense posture. 

If our ships could have the option of burning the most 
readily available fuel, this reduction in readiness can be 
avoided by switching to coal, a fuel source in great abundance 
trietne U.S. 

Paramount in the engineering of a successful multi-fuels 
propulsion plant is the development of a marinized fluid bed 
combustor. The design of heat exchangers and other turbo- 
machinery used in tne combined cycle are well within the state 
Oi the art of present technology. 

It is thus that the author concludes the following: 

1) The design of a modern coal burning combined cycle is 
Within the reach of present technology. 

2) This combined cycle would offer high efficiency at 
both rated load and part load due to the unique controlling 
aspect of the Closed Brayton cycle. 

3) Weight and volume of a Closed Brayton combined cycle 
would be greater than that of an equivalently rated COGAS type 


system. 


4) Coal offers the benefit of a readily available, 


> 





relatively cheap fuel at the expense of increased ship size 
and acquisition cost due to lower BTU heat value per ton when 


compared to oil. 


The author makes the following recommendations: 

1) Development of a land based test site of a Closed 
Brayton combined cycle equipped with a fluid bed heater should 
be conSidered. 

2) A pressurized fluidized bed combustor suitable to 
marine applications should be developed. 

3) A shipboard prototype system should be integrated 
ona U.S. Navy ship to attempt to quantitatively define the 


shipboard impact of the system. 
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APPENDIX. A 


The computer program performs a thermodynamic analysis 


and component volume computation for various specified 


Brayton cycle compressor pressure ratios. 


The program accepts aS input the following quantities: 


Vc (P) 
/+ (P) 


ER 


Deis) 


Bey = 


(= 
{ 


 WHB 


SAT 


Power 
ie 
Fag 


ic Ph 


/* (is) 


polytropic compressor efficiency (Brayton 
cycle) 


polytropic turbine efficiency (Brayton 
cycle) 


regenerator effectiveness 


isentropic pump efficiency (Rankine 
cycle) 


total pressure drop around Brayton cycle 
loop 


efficiency of the waste heat boiler 


compressor inlet temp. (Brayton cycle) 
turbine inlet temp (Brayton cycle) 
steam condensing temp. 


min. temp diff. between streams in WHB 
(at ies Fey alte! a 


min. Saturation pressure in the Rankine cycle 


net power out of the combined cycle 


compressor inlet temp. (heater cycle) 
turbine exit temp. (heater cycle) 


polytropic compressor efficiency (heater 
cycle) 


isentropic turbine efficiency (Rankine 
eycle) 
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Ye(P) i - polytropic turbine efficiency (heater cycle) 
g - equivalence ratio 


WP - working pressure in the Brayton cycle. 


Output is as shown in Appendix B. 

The main program performs the thermodynamic analysis of 
the closed Brayton cycle and the Rankine cycle. Functions per- 
formed by the various subfunctions are as follows: 

1) Subroutine HEATER performs a thermodynamic analysis of 


the heater cycle. 





2) Subroutine DHD computes = for the regenerator, re- 
cuperator, cooler and heater. 
3) Subroutine DPP calculates for the regenerator, re- 
Suwerauor, cooler, heater and boiler. 
4) Subroutine TURVOL finds turbine volumes. 
5) Subroutine COMPVOL finds compressor volumes. 
6) Subroutine REGEN computes the regenerator volume. 
7) Subroutine RECUP computes the recuperator volume. 
8) Subroutine COOLER computes the cooler volume. 
9) Subroutine COMBUSTOR computes the heater volume. 
10) Subroutine WHBVOL computes the waste heat boiler volume. 
11) Subroutine CONDVOL computes the condenser volume. 
12) Subroutine STMTAB is a computerized version of the steam 


tables. Reproduced from ref. 15. 
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10 
15 
20 


30 
40 


50 
60 


DIMENSION CR(30) ,SATLIQ(3) ,STEAM(3) , TPP(12) ,WR(11) 
REAL NCP NTP NPIS,NTIS,NWHB,NT,L1,L2 ,KO 
REAL NCPH,NTPH,MB,MS,MA,MF NHN, IT,KIL,NPP 
REAL LHV,M,M1,M2,NTU 
COMMON /STM/PRESS,TEMP ,QUAL,SATLIQ, STEAM 
COMMON /HTR/NH,PRH,DPPH,T2H, T3H, T4H, T5H, TOH MA, MF 
COMMON /DIA/DHDREG,DHDCLR,DHDHTR ,DHDREC 
COMMON /PRES/DPPWHB ,DPPREG ,DPPCLR ,DPPREC ,DPH1 ,DPH2 ,DHDHT 
COMMON /COMB/DPPRG,DPPRC ,DPPHE ,DHDHR ,DPHE2 
READ (5,10) NCP,NTP,ER,NPIS,NTIS,DP,NWHB,T1,14,T1S,DT,PMIN,K 
FORMAT (7F4.3,5F7.1,12) 
READ(5,15) PWR,TLH,T7H,NCPH,NTPH,PHI, WP 
FORMAT (F8.1,2F6.1,3F5.3,F6.1) 
READ(5,20) (CR(1),I=1,K) 
FORMAT (<K>F6.2) 
R= 068544 
mee z00 1=1,k 
PR=CR(I) 
Z=CP1(T1)+(R/NCP) *LOG(PR) 
T2=EXP(Z/ 2475) 
BALLOITER (T2,Z) 
TR=1/(1-DP)/PR 
Z=CP 1( T4)+R*NTP*LOG( TR) 
T5=EXP(Z/.2475) 
CALL ITER (T5,Z) 
etme .Gl. 15) GO TO 30 
T6=T5-ER*(T5-T2) 
T3=T2+ER* (T5-T2) 
GO TO 40 
T6=T5 
T3=T2 
Li=Tl 
e7-i22.0 
fiete2 .Gi. 16) GO 10 50 
GO TO 60 
L2=16-(16-T1)/10. 
TPP(1)=L2 
WR1=0. 
t=(e2-1 1) /10. 
DO 70 J=1,11 
TSF=TPP(J)-DT 
ieS=10-0 1 
T=TSF-460 
CALL STMTAB (2,T,0.) 
PSF=PRESS 
HSF=SATLIQ(2) 
T=T3S-460 
CALL STMTAB (7,PSF,T) 
H3S=STEAM(2 ) 
S3S=STEAM(3) 
H=111S=460 
61 





68 
70 
is) 


76 


uy 


CALL STMTAB (2,T,0.) 
P4S=PRESS 
V1S=SATLIQ(1) 
HIS=SATLIO(2) 
S1S=SATL10Q(3) 
CALL STMTAB (9,P4S,S3S) 
H4SIS=STEAM(2) 
H4S=H3S-NTIS*(H3S-H4S IS) 
CALL STMTAB (8,P4S,H4S) 
X=QUAL 
V4S=STEAM(1) 
H2SIS=H1S+V1S*(PSF-P4S ) 
H2S=H1S+(H2SIS-H1S) /NPIS 
T2S=T1S+15. 
B=NWHB* (CP2(1T6)-CP2(TPP(J)))/(H3S-HSF ) 
Z=ACP(TPP(J))*TPP(J )-B* (HSF-H2S ) /NWHB 
T7=TPP(J) 
CALL ITER2 (17,Z) 
WR(J)=B* ( H3S+H1S-H4S-H2S ) 
IF(J .£Q. 1) GO TO 68 
tei iieelh. 11) GO TO 75 
IF(PSF .LT. PMIN) GO TO 75 
IF(WR(J) .GT. WR1) GO TO 68 
me eeGhs <07) GO 10 75 
WRL=WR(J) 
TPP(J+1)=TPP(J)-C 
CONTINUE 
GO TO 76 
IF(ABS(TPP(J)-TPP(J-1)) .LT. .1) GO TO 76 
L1=TPP(J) 
PZ=PPtd-1) 
GO TO 60 
WB=CP 2( T4)+CP2(T1)-CP2( 75) -CP2( 72) 
QIN=CP2( T4)-CP2(T3) 
NT=(WB+WR(J) )/QIN 
MB=.707*PWR/(WB+WR(J) ) 
MS=B*MB 
QB=MB*QIN 
CALL HEATER (QB,T1H,T7H,NCPH,NTPH,PHI, 13,14) 
CALL DHD(PR) 
DHDHT=DHDHTR 
er 
CALL DPP(T6,T7,PR,DP,ER,T2,7T3,DPPH,T4H,15H,1T4,PRH,WP,MB,MA 
MF ,K) 
G=1.33 
N=G/(G-(G-1.)*NTP) 
P=WP* (PR-DPPREG) 
V=1./ARHO(T4,P) 
VMNTURB=0.0 
CALL TURBYOL (N,MB,V,TR,1T4,VMNTURB) 
G=1.39 
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N=G/(G-(G-1.)*NTPH) 

P=14,7* (PRH-DPPREC) 

V=1. /ARHO(T6H,P) 

M=MA+MF 

VHTRTURB=0 .0 

TRH=1./(1.-DPPH)/PRH 

CALL TURBVOL (N,M,V,TRH,T6H,VHTRTURB) 

T=T1S/T3S 

P=P4S/PSF 

N=1./(1.-L0G(T)/LOG(P) ) 

VSTMTURB=0.0 

CALL TURBVOL (N,MS,V4S,P,T3S,VSTMTURB) 

G=1.4 

N=G/(G-(G-1.)/NCP) 

V=1./ARHO(T1, WP) 

VMNCOMP=0.0 

CALL COMPVOL (N,MB,V,PR,VMNCOMP ) 

V=1. /ARHO(535. ,14.7) 

VHTRCOMP=0.0 

CALL COMPVOL (N,MA,V,PRH,VHTRCOMP ) 

D=.125 

VREGEN=0.0 

CALL REGEN (T2,T3,T5,T6,WP,D,MB,PR,DPPREG,DHDREG, VREGEN) 
VRECUP=0.0 

CALL RECUP (T2H,T3H,T5H,T6H,PRH,DPPH,D,MA,MF .DPPREC ,VRECUP) 
D=.25 

NeER=0.0 

CALL COOLER (17 ,T1,WP,D,NTP,MB,DHDCLR ,DPPCLR,VCLR) 
DPPRG=DPPREG 

DPPRC=DPPREC 

DPPHE=DPPH 

DHDHR=DHDHTR 

DPHE2=DPH2 

VHTR=0.0 

CALL COMBUSTOR (T3,T4,T4H,T5H,WP,PR,D,MB,PRH,MA,MF ,VHTR) 
VWHB=0.0 

CALL WHBVOL (MB,MS,T6,T7 ,NWHB,H2S ,H3S ,DPPWHB ,VWHB) 
VCOND=0.0 

CALL CONDVOL (MS,X,VCOND) 
VOLT=VMNTURB+VHTRTURB+VS TMTURB+VMNC OMP+VHTRCOMP+VREGEN+VRECUP 
+VCLREVHTR+VWHB+VCOND 

IF(K .£EQ. 1) GO TO 78 

IF(VOLT .GT. VOL1) GO TO 81 

IF((DPPCLR-.002) .LT. 0.0) GO TO 81 

K=K+1 

VOL1L=VOLT 

V1=VCLR 

V2=VWHB 

GO TO 77 

VOLT=VOL1 


VCLR=V1 
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VWHB=V2 
WRITE(6,82) PR,WP,DP,MB 
82 FORMAT('1', ‘BRAYTON CYCLE: PRESSURE RATIO=',F5.2/1X, ‘WORKING 
* PRESSURE=',F5.1/1X,'TOTAL PRESSURE DROP=',F3.2/1X, 'MASS FLOW 
* (1bm/s)=',F6.2//10X, ‘STATE TEMP(deg R)') 
WRITE(6,83) Tl 
83 FORMAT (12X,‘'1',6X,F6.1) 
WRITE(6,84) T2 
84 FORMAT (12X,'2' ,6X,F6.1) 
WRITE(6,85) T3 
85 FORMAT (12X,'3',6X,F6.1) 
WRITE(6,86) 14 
86 FORMAT (12X,'4' ,6X,F6.1) 
WRITE(6,87) T5 
87 FORMAT (12X,'5' ,6X,F6.1) 
WRITE(6,88) T6 
88 FORMAT (12X,'6',6X,F6.1) 
WRITE(6,89) TPP(J) 
89 FORMAT (12X, 'PP' ,5X,F6.1) 
WRITE(6,90) T7 
90 FORMAT(12X,'7',6X,F6.1//5X,'COMPONENT DPP DHD VOL(cu ft)') 
WRITE (6,91) DPPREG,DHDREG,VREGEN,ER 
91 FORMAT(7X, 'REGEN' ,4X,F4.3,1X,F4.3,1X,F7.1,2X, 'EFF=',F4.2) 
WRITE (6,92) VMNTURB ,NTP 
92 FORMAT(6X, 'TURBINE' ,15X,F5.1,2X, 'EFF(P)='F4.2) 
WRITE(6,93) DPPCLR,DHDCLR,VCLR 
93 FORMAT(6X, 'COOLER' ,4X,F4.3,1X,F4.3,3X,F5.1) 
WRITE (6,94) VMNCOMP NCP | 
94 FORMAT(4X, ‘COMPRESSOR’ ,14X,F5.1,2X, 'EFF(P)=',F4.2) 
WRITE (6,95) DPH1 
95 FORMAT(6X, 'HEATER' ,4X,F5.4) 
WRITE (6,96) DPPWHB ,NWHB 
96 FORMAT(6X, 'BOILER' ,4X,F4.3,15X, 'EFF=',F4.2/) 
WRITE(6,100) PRH,DPPH,MA,MF 
100 FORMAT(1X, ‘HEATER CYCLE: PRESSURE RATIO=',F5.2/1X, ‘TOTAL 
* PRESSURE DROP=',F3.2/1X,'AIR FLOW(1bm/s)=',F6.2/1X, ‘FUEL 
* FLOW(1bm/s)=',F6.2//10X, ‘STATE TEMP(deg R)') 
WRITE(6,101) T1H 
101 FORMAT(12X, '1h' ,5X,F6.1) 
WRITE(6,102) T2H 
102 FORMAT(12X, '2h' ,5X,F6.1) 
WRITE(6,103) T3H 
103 FORMAT(12X,'3h' ,5X,F6.1) 
WRITE (6,104) T4H 
104 FORMAT(12X, '4h' ,5X,F6.1) 
WRITE(6,105) T5H 
105 FORMAT(12X, '5h' ,5X,F6.1) 
WRITE(6,106) T6H 
106 FORMAT(12X, ‘6h' ,5X,F6.1) 
WRITE(6,107) T7H 
107 FORMAT(12X,‘'7h' ,5X,F6.1//5X, ‘COMPONENT DPP DHD VOL(cu ft) ' 
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110 
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Piz 
3 
20 


Lat 
Bae 
WAS 
124 
2s 
128 
129 
130 
135 
136 


Se 


1338 
200 


ERH=(T3H-T2H) /(T5H-T2H ) 


WRITE(6,110) DPPREC ,DHDREC,VRECUP,ERH 

FORMAT(7X, 'RECUP' ,4X,F4.3,1X,F4.3,1X,F7.1,2X, 'EFF=',F4.2) 
WRITE (6,111) VHTRTURB ,NTPH 

FORMAT(6X, 'TURBINE' ,15X,F5.1,2X, 'EFF(P)=',F4.2) 
WRITE(6,112) DPH2 ,DHDHTR,VHTR 

FORMAT(6X, 'HEATER' ,4X,F4.3,1X,F4.3,3X,F5.1) 

WRITE (6,113) VHTRCOMP ,NCPH 

FORMAT(4X, 'COMPRESSOR' ,14X,F5.1,2X, 'EFF(P)=',F4.2/) 
WRITE (6,120) PSF ,MS,B 

FORMAT(1X, 'RANKINE CYCLE: SATURATION PRESSURE=',F6.1/1X, 


* "STEAM FLOW(1bm/s)=' ,F6.2/1X, ‘FLOW RATIO=',F6.4//10X, 'STATE 
* TEMP(deg R) PRESSURE(psi)') 


WRITE(6,121) T1S,P4S 
penile x. 0S ,59%,F0,),8X,F6.1) 

WRITE(6,122) T2S,PSF 
BORMAN(IZX, «25 ,5X,FO.1,8X,F6.1) 

WRITE(6,123) TSF ,PSF 

BORMAN I2X, SF: ,5X,FO.1,8X,F6.1} 

WRITE(6,124) T3S,PSF 

FORMAT(12X,'3s' ,5X,F6.1,8X,r6.1) 

WRITE(6,125) T1S,P4S,X 

FORMAT(12X, '4s' ,5X,F6.1,8X,F6.1,2X, ‘QUAL=' ,F5.3/) 
WRITE(6,128) VWHB 

FORMAT(5X, ‘COMPONENT VOL(cu ft) '/6X, ‘BOILER’ ,5X,F6.1) 
WRITE(6,129) VSTMTURB ,NTIS 

FORMAT(6X, 'TURBINE' ,5X,F5.1,2X, ‘EFF(IS)=' ,F4.2) 
WRITE(6,130) VCOND 

FORMAT(5X, ‘CONDENSER’ ,4X,F5.1/) 

Mete (6,135) NT 

FORMAT('1L' , "SUMMARY: THERMAL EFF=' ,F5.3) 
WRITE(6,136) NH,PHI 

FORMAT(9X, 'HEATER EFF=',F5.3,' (PHI=',F4.2,')') 
SF C=MF*3600./PWR 

WRITE(6,137) SFC,VOLT,PWR 

FORMAT (9X, 'SFC(1bm/HP-hr)=' ,F5.3/9X, ‘VOL TOTAL(cu ft)=',F6.1/ 


*9X, ‘POWER(HP)=' ,F7.1) 


RANKINE=WR (J) /WB 

WRITE(6,138) RANKINE 

FORMAT (9X,'WORK FRAC BY STM CYCLE=' ,F5.3) 
CONTINUE 


= SuOP 
END 
FUNCTION CP1(T) 
CP1=.2475*L0G(T)-(3.759E-5*T )+(2.553E-8*T**2 )-(4.103E-12*T*%3 ) 
RETURN 
END 
FUNCTION CP2(T) 
CP2= ~24/5*T- ( 1 ~879E-5* T**2 )+( l P 702E-8* T**3 ) -(3 : O78E-12* 7**4 ) 
RETURN 
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10 


20 


END 

FUNCTION AMU(T) 

Hi 5 1/9 

AMU=TK**1 .5*1.E-6/ (TK+110.4) 

RETURN 

END 

FUNCTION ARHO(T,P) 

Rao. 30 

ARHO=P*144./R/T 

RETURN 

END 

FUNCTION AK(T,P) 

A= .6325E-5 

B=245.4 

C=12. 

i= 3.7 1/9. 

Pee ATK**.5) /{1 .+B*10.**(=1.*C/1K)} /1K) *418.6 
AK=( AKO+(1.464E-5* (16 .02*ARHO(T,P) )**1.23))*.578/3600. 
RETURN 

END 

FUNCTION APR(T,P) 

CP=ACP(T) 

XMU=AMU (T) 

XK=AK(T,P) 

APR=CP *XMU/XK 

RETURN 

END 

FUNCTION ASIGMA(T,P) 

ASIGMA=17.214* (6228. *AK(T,P) )**(-1.41) *(1.46*AMU(T) )**1.23% 


*(1./APR(T,P)**.564)*(1./(16 .0O2*ARHO(T,P))**.41) 


RETURN 

END 

FUNCTION ACP(T) 

Mibeme co) J=3, | 9E~o*1+5.1L06E-8*1**2-1.231E-11*1**3 
RETURN 

END 


SUBROUTINE ITER (T,X) 

T1sT 

T=EXP( (X+.2475*LOG(T)-CP1(T))/.2475) 
IF(ABS(T-T1) .LT. 1.E-2) GO TO 20 
GO TO 10 

RETURN 

END 

SUBROUTINE ITER2(T,X) 

T1=T 

T=X/ACP(T) 

IF(ABS(T-T1) .LT. 1.E£-2) GO TO 20 
GO TO 10 

RETURN 

END 
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20 


SUBROUTINE HEATER (QB,T1H,T7H,NCPH,NTPH ,PHI,13,T4) 
REAL LHV,NH,MF MA ,NCPH,NTPH,N 

COMMON /HTR/NH,PRH,DPPH, 12H, 13H, T4H,T5H, TOHMA MF 
LHV=12900. 

FA=.0922 

H7H=ACP (17H )*T7H 

H1H=ACP(T1H)*T1H 

NH=1-( (1+FA*PHI) *H7H-H1H ) / (FA*PHI*LHV) 
MF =QB/ (NH*LHV ) 

MA=MF / (FA*PHI) 

T2H=0 .0 

T6H=0.0 

CPT=ACP( ( T6H+T7H) /2. ) 

CPC=ACP( (TLH+T2H)/2. ) 

N=(1+FA*PHI )*CPT*NTPH* (CPT*T7H-CPC*T1H ) 
D=CPC*T1LH* (CPC /NCPH-CP T*NTPH) 
PRH=(N/D)** (NCPH/CPC ) 
TH2=TLH*PRH** (CPC /NCPH) 

N=(1+FA*PHI) *CPT*T7H+CPC*T1H* (PRH**(CPC/NCPH )-1) 
D=(1+FA*PHI)*CPT*T7H*PRH** (CPT/NTPH) 
DPPH=1-(N/D)**(1./(CPT*NTPH) ) 
TRH=(1-DPPH)*PRH 
TH6=T7H* (TRH) ** ( CPT*NTPH ) 
IF(ABS(T2H-TH2) .LT. .1) GO TO 20 
IF(ABS(T6H-TH6) .LT. .1) GO TO 20 
T2H=TH2 

T6H=TH6 

GO TO 10 

T2H=TH2 

T6H=TH6 

T5H=T3+100. 
T4H=T5H+QB / (ACP (T5H /2+1300. )*(MA+MF) ) 
T3H=T2H+(1+MF /MA)* ( TSH-T6H) 

RETURN 

END 

SUBROUTINE DHD(PR) 

COMMON /DIA/DHDREG,DHDCLR ,DHDHTR ,DHDREC 
DHDREG= .557+.0933*PR 

IF(DHDREG .GT. .95) DHDREG=.95 
DHDCLR=.076+.092*PR 

DHDHTR= .214+.048*PR 

DHDREC=.65 

RETURN 

END 


SUBROUTINE DPP(T6,T7,PR,DP,ER,T2,13,DPPH,TH4,THS ,T4 ,PRH,P,M,MA 


eer ,K) 


REAL M,MA,MF M1, M2 
COMMON /PRES/DPPWHB ,DPPREG,DPPCLR ,DPPREC ,DPH1, DPH2 ,DHDHT 
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DTG=T6-T/ 

DP PWHB=FLOAT (K)* .002 
DPPREG=(.135-.0284*PR)*DP/.1 
IF(DPPREG .LT. .01) DPPREG=.01 
mize. £0. WS) DPPREG=0.0 
DPH2=( .0108*PR-.01)*DP/.1 
f=(13ti4)*.5 
TT=(TH4+TH5)* .5 

P1=P*(PR-DP) 
P2=14.7*(PRH-DPPH) 

Ml=M 

M2=MA+MF 

XMU1=AMU(T1) 

XMU2=AMU(TT ) 
XRHO1=ARHO(T1,P1) 
XRHO2=ARHO(TT ,P2) 

B=(M1 /M2)**(-1. 75) *(XMU1L/XMU2 ) ** (- .25) *(XRHOL/XRHO2 ) *(P1/P2) 
DPH1=DP H2*DHDHT**3/B 
DPPCLR=D P-DPPWHB-DPPREG-DPH1 
DPPREC=DPPH-DPH2-.1 

RETURN 

END 


SUBROUTINE TURBVOL (G,M,V,TR,T,VOL) 
REAL M,N,KI1,IT 
SF=4.0 
UT=850. 
PHIT=0.6 
VA=0.9 
CL=2.0 
PI=3.14159 
RHOM=15.12 
A=9 .63E6 
B=-95.14 
ALPHA=1.2 
Al=2.5 
n= .307 
K11=3.39£-3 
THETA=0.8 
RT=SORT(M*V/ (PHIT*UT*P I*(1.-VA**2) ) ) 
VB=SQRT(1.-TR**(1./G)*(1.-VA**2) ) 
C1=0.2* (AL+1. )*1. /V¥(1. -VA**2) SUT**2*RT**2*X 1/KIL 
C2=2.*(A-B*T) /(SF*ALPHA) 
C3=2.*B*T*(1.-VA**2 ) **(G-1. ) /(SF*ALPHA) 
C4=-1.*THETA*RHOM*UT**2 /2, 
A= VA-VB ) /20. 
SUM=0.0 
E00 
Y1=VA 
DO 100 J=1,21 
N=C1*((1.-Y1)/(1.+Y1)) *PHIT**2+Y1**2 
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50 
100 


Daeanesanea Yam 2)o-| eeG)+C4* (1. = 1**2 ) 
Pen ee Veer Gleecey te -4 PHID<*2+Y 1**2)/(1.4+Y1)=(1.+Y1} 


pone 2 = (1. -¥1)*(PHI1**2+Y¥ 1**2) ) 


BO eels (to-y 1**2)**(-1.*G)*(G-1. )*C3-2. *Y1*C4 
X2=.5*(1./(N/D)**.5)*( (DN*D-DD*N ) /D**2 ) 
Bein 20. 0.0) GO 70 50 
SUM=SUM+(X1+X2) /2. 
X1=X2 
Y1=Y1-H 
CONTINUE 
TT=H*SUM 
wOleec. ~P 1s Ri<*2*CL*1T 
RETURN 
END 


SUBROUTINE COMPVOL (G,M,V,PR,VOL) 
REAL M,N,KIL,IT 

SF=4.0 

UT=850. 

PHIT=0.5 

VA=0 .6 

l=2.0 

PI=3.14159 

RHOM=15.12 

ALPHA=1 .2 

pi=12.0 

X1=.05 

K11=4,936-5 

THETA=0.8 

SIGMAY=12.54£6 

RT=SORT (M*V/ (PHIT*UT*PI*(1.-VA**2) ) ) 
VB=SQRT(1.-PR**(-1./G)*(1.-VA**2) ) 
C1=.2*(AL+1.) *1./V* (1. -VA®*2) XUT* *2*RT**2*PHIT**2*X1 /KI1 
C2=-1* THETA*RHOM*UT**2 /2. 

C3=2.*SI GMAY/(SF*ALPHA) 
H=(VB-VA)/20. 

SUM=0 .0 

X%1=0.0 

YL=VA 

meo0| J=1,21 
N=C1*((L.-Y1)/(1L.+Y1) )*(PHIT*#*2+.49*Y 1**2) / (PHIT**2+.7225*Y 1**2 


rm) 


p-es+C2*(1.-¥1**2) 
ene neh set. 49° 1 <2 )/(PHIT**2+./225*1 1**2 )*(1.4Y1)-(1 


eet lie (= 2)" (PHI ix =2>).4 9%) 1**2) / (PHIT**2+.7225*Y 1**2 } + 
meson 1 (1 Yi) /(1.+Y1) /(PHIT**2+.7225*Y 1**2 )-1 .445*Y1*(1.-¥1) / 
See WAP Hi) **2h 4 9*Y 1**2)* (PHI T**24 .7 225*Y 1**2)**(-2) ) 


DD=-2. *C2*Y1 

X2=.5*(1./(N/D)**.5)*( (DN*D-DD*N) /D**2 ) 
Maixl -E0. 0.0) GO TO 50 

SUM=SUM+ (X1+X2)/2. 
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99 


X1=X2 

Y1=Y1+H 

CONTINUE 

IT=-H*SUM 
wOle=2. =P 1*RE**2*CL*1T 
RETURN 

END 


SUBROUTINE REGEN (12,13,15,16,P,D0,MB,PR ,DPPREG,DHDREG, VOL) 
REAL MB 

T=715/2+1T6/2 

Pi=P*PR 

SIGMA=ASIGMA(T,P1) 

ieii2 .cQ. 13) GO TO 99 

DTLM=T3-T2 

Gevite=1./ (6003. *P1)**.41*{ .0254*D)**1 .18/ (.56*DTLM)**1 .41 
DH1=CP2(T5)-CP2(T6) 

PT=3.14159 

AREG=P 1/4 .*SIGMA*GAMMA* (MB* .454) *(DH1*2 326. )**1.41* (1.+DHDREG) 


ee) * (1 tPRAX( +2) *DEDREG** (-3) )** 41735 .3 


VOL=AREG*DPPREG** (-.41) 
RETURN 
END 


SUBROUTINE RECUP (12H,T3H,T5H,T6H,PRH,DPPH,D,MA,MF ,DPPREC ,VOL) 
REAL MA,MF 

eto /2.+16H/2. 

PC=1./(1.-DPPH) 

P=14./7*PRH 

SI GMA=ASIGMA(T,P) 

DTLM=T3H-T2H 

GAMMA=(6803. *P) **(~.41) *(.0254*D )**1 .18*( .56*DTLM) **(-1.41) 
DH1=CP2( TSH )-CP2( T6H) 

PI=3.14159 

DHDREC= .65 

CM=MA+MF 

AREC=P 1/4. *SIGMA*GAMMA* (CM* £454) *(DH1*2326. )**1.41* (1 .+DHDREC) 


#kKO AL* (1, +DHDREC** (-3) *PC**2)** 41435 .3 


VOL=AREC*DPPREC**(- .41) 
RETURN 
END 


SUBROUTINE COOLER (17,T1,P,D,NTP,MB,DHDCLR,DPPCLR,VOL) 

REAL NPP ,NTP,MB 

IF(T7 .LT. T1) GO TO 99 

TA=17/2.+T1/2. 

SIGMA=ASIGMA(TA,P ) 

TC1=530. 

TC2=540. 

DTLM=(T1-TC1-T7+TC2) /LOG( (T1-TC1) /(17-TC2) ) 

GAMMA= (6803. *P)** (-.41)*( .0254*D)**1.18*( .S6*DTLM)** (-1.41) 
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DH1=CP2(T7)-CP2(T1) 

P1=3.14159 

TW=TC1/2.+7C2/2. 

RK=AK(TA,P)/1.£-4 

RMU=AMU(TA) /6.4E-4 

RM=(T7-T1) *ACP(TA)/(1.*10. ) 

RPR=APR(TA,P) /6.4 

NPP=.80 

RRHO=ARHO(TA,P) /64 .0 

RSV=ARHO(T7 ,P) /64.0 

ALPHA=RK*RMU** (-.8)*RM**, 8*RPR** 4 

CB=RM**(-2.75) *RMU**(-.25) *RRHO*RSV/NTP/NPP 
ACLR=P1/4.*SIGMA*GAMMA* (MB* .45 4) * (DH1*2326. )**1 .41* (1.+DHDCLR) 
** (1. +ALPHA*DHDCLR) **1 .41* (1 .+C B*DHDCLR**(-3) )**.41*35.3 
VOL=ACLR*DP PCLR** (-.41) 

RETURN 

END 


SUBROUTINE COMBUSTOR (173,T4,T4H,T5H,WP,PR,D,MB,PRH,MA,MF VOL) 
REAL MB,MA,MF ,M 
COMMON /COMB/DPPRG,DPPRC ,DPPHE ,DHDHR ,DPHE2 
T=T3/2+7T4/2 
P=WP* (PR-DPPRG) 
SIGMA=ASIGMA(T,P) 
DTLM=(T5H-T3 -T4H+T4 ) /LOG( (T5H-T3) /(T4H-T4) ) 
GAMMA=(6803.*P)**(-.41)*( .0254*D) **1 .18*( .56*DTLM)**(-1.41) 
DH1=CP 2(T4)-CP2(T3) 
PI=3.14159 
TH=T4H /2+T5H /2 
TRH=(1.-DPPHE)*PRH 
PH=14.7* (PRH-DPPRC) 
RK=AK(T,P) /AK(TH,PH) 
RMU=AMU(T) /AMU (TH) 
M=MA+MF 
RM=MB/M 
RPR=APR(T,P)/APR(TH,PH) 
ALPHA=RK*RMU**(- .8) *RM**, 8*RPR** 4 
AHTR=P1/4.*S1GMA*GAMMA* (MB* .454) * (DH1*2326.)**1 .41*(1.+DHDHR) * 
*  (1.+ALPHA*DHDHR) **1 .41*DHDHR**(-1.23) *35.3 
VOL=AHTR*DPHE2** (-.41) 
RETURN 
END 


SUBROUTINE WHBVOL (MB,MS,T6,T7 ,NWHB,H2S ,H3S ,DPPWHB ,VOL) 
REAL MB,MS,NWHB ,NTU 

fievo/2t1//2 

CP=ACP(T) 

CPS=2.0 

CR=MS*CPS/(MBX*CP ) 

PetcR .GT. 1.) GO TO 10 

CMIN=MS*CPS 





10 
20 


GO TO 20 

Grae /CR 

CMIN=CP*MB 

PS=4.0 

ER=(NWHB*CR-1.)/(NWHB-1. ) 
tee ERUUNL. /PSTYACR ER (L./ 
p 

DH=H3S-H2S 

VOL= .00045*CMI N*NTU*MS*DH*DPPWHB** (- .41) 
RETURN 

END 


SUBROUTINE CONDVOL (MS,X,VOL) 
REAL MS 

VOL=MS*X*27.3 

RETURN 

END 
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SUBROUTINE STMTAB (NI,XI,TI) 


DIMENSION SATLIQ(3) ,STEAM(3) ,COEFFT (237) 
COMMON /STM/PRESS, TEMP ,QUAL,SATLIQ, STEAM 


pera (COEFF I(1) ,1=1,30) 
modo Ss7lE+O3, .60464261E+01, . 
- .19513682E+00 ,- .52279952E+00, . 
- 50625849E+00, .72595/720E+00, 
.55592942E+03, .20000000E+01,-. 
.26548914E-04, .91349306E-04, 
Pou OGII9IETO0, .36241343Et01, 
-45000000E+03, .55000000E+03,-. 
pare (COEFFI(1),l= 31, 60) / 
- .32893178E-14, -.11734466E-11, 
molosoocIE-O2, .69519265E+01, 
~.64272696E-21, .12/704366E-18,-. 
~ .64251564E-10, .44797322E-07,-. 
.33953897E+01, .43010158E+03, 
9275654 3E-17 ,- .95470656E-15,-. 
.39936598E-08,-.99567313E-06,-. 
ATA (COEFFT(1),1= 61, 90) ih 
pred) 22907E—21, .66553573E-19, 
Wa 4555909E-12,— .18335302E-10,-. 
.98512754E+00, .25000000E+03, 
m55/ 1S0E-22, .33225506E-19, 
-./0679233E-14, .90114201E-11, 
.30734203E-04, .21762691E-01, 
- .33608032E-23, .85819792E-20, 
pia (COEFFI(I) ,I= 91,120) i 
.15587770E-07, .81264440E-05, 
55000000E+03, .94743147E-18, 
.44979672E-12,-.51337380E-10, 
»12399614E-02, .12798551E+01, 
~12945371E-21, -.65873139E-19,-. 
»17548001E-12, .23202424E-09, 
-10150974E+01, .21848036E+03, 
ATA (COEFFT(1) ,1=121,150) /-. 
.69217855E-16,-.18345822E-14,-. 
59598410E-08, .45575915E-06, 
~25000000E+03, .50084539E-23,-. 
-49239634E-17, .16250351E-13,-. 
-.85546496E-06, .14285709E-02, 
-L9000000E+01,-.11085077E+05,-. 
-.18796109E+03, .15699797E+03,-. 
ATA (COEFFT(1),1=151 ,180) / 
-11199490E+04, .15000000E+01,-. 
-62379937E+05, .29474498E+05,-. 
.35122466E+04, -.12340833E+04, -. 
-30000000E+01, .13000000E+04, 
.52122736€+03,-.11682947E-01,-. 
.12340954E-04, .13555136E-03,-. 
.18000000E+01, .23822091E+02, . 


+ + + + + 4H 
Ow «wv ww ww ww © w Ce © © we © © © 


+ + t+ + t 4 


CJewe «w we we we & w& 


+ + + + * 4 F 


70000000E+01, 
91148958E-01, 


.14585694E+02, 


74133440E-06, 


.15679922E-02, 
.44960616E+02, 


98/36888E-22, 


.30041189E-17, 
ead Ge 507. 
.25000000E+03, 


96976163E-16, 
28104756E-04, 


poe ol 1ELOS., 


29004677E-12, 
Sg) 792E-06. 


55000000E+03, 
.47725605E-17, 


31160333E-06, 


47000000 +03, 
3007557 5E-18, 
58149642F-09, 
.25000000E+03, 
.25729365E-18, 
.60506309E-13, 
.17002639E-01, 
.24971478E-16, 
.49178986E-08, 
5492737 7E+03, 


13817574E-17, 


.36145967E-06, 
-45000000E+03, 


21134438E-20, 
66417458E-12, 


.12325099E-02, 


15659070E-21, 
74068320E-12, 


.36753405E+00, 


15994587E+05, 
20213342E+03, 


.15040455E+03, 


14863585E+06, 
18907586E+05, 
68928827E+02, 


,70848990E+01, 


12001588E+00, 
1509931 9E-03, 
18776649E+03, 


Sind Oh’ 


= 


Vmeccuecee \E+04.. .65000000Et0S 
.18552680E+00 
.61843682E-01 
.12404319E+03 
1332482 3E-06 
-25032530E-01 
.17926732E+03 
.19301230E-19/ 
.42113669E-15 
./8419653E-05 
.20575548E-23 
.83370000E-13 
.17490541E-01 
.18444779E-19 
~35353292E-10 
~/8143274E+00/ 
2184803 6E+03 
.53349163E-14 
.11956354E-03 
~55000000E+03 
./3898651E-15 
.86712981£-07 
-7/8459129E-25 
~31591904E-15/ 
.13684599E-10 
-45000000E+03 
3830258 8E-14 
~48326895E-05 
25000000E+0 3 
.62520106E-14 
.12506461E-03 
~55000000E+03/ 
.15327038E-18 
.14589563E-10 
./4969797E+00 
.46099119E-18 
.10479899E-08 
.1 7000000E+01 
.91407688E+04 
-48650176E+02/ 
.19050547E+03 
.13414819E+06 
.5/85887 6E+03 
.12036836E+04 
.16507788E+02 
.20411991E+01 
3000000 0E+01 
.33080504E+03/ 





1000 
1010 


1020 
1030 
1040 
1050 
1060 
1070 


1080 


1090 


1100 
1110 
1120 
MING 


1140 
eer’ 


te + OF 


Cw ww ~~ ww bod bod ww 


+ + + + 


PW AMCGOEFF IV 1 )41=181, 210) /-.32141408E+01, .42449101E+03 
.16400835E+04, .18746657E+02, .63978878E+02, .20719841E+04 
~13000000E+04, .18000000E+01,-.64369800E-05, .96418803E-02 
.38879248E+01, .11034060E-05, .47225001£-03,-.94113600E+01 
po oe logy SE-05,— c0G92120E-02, .10951922E+01, .13000000E+04 
mrovOUOUOETOL,—_3s3101206E-03, .Z21661916E+01, .53500Z291E+03 
.77236618E-03, .36753650E+00,-.25450121E+03, .25249707E-02 
~.16427106E+01, .61477766E+03, One 0.0 
pany (COEFFT( 1) ,1=211, 237) / Cameo. COO Ore sO 0 
Teno U 0. 0.0,0.0,0.0,0.0,0.0,0.0,0.0, .J0000000E+03 
L4000000E+01,-.15478259E-04, .19351556E-01, .24698909E+01 
-23329644E-04, .11493791E-02,- .13139208E+02, .58860219E-05 
Smee go2S99E-Ul, .solOLIGUE+01/ 


DELTA1=0 .0 

ERRORL=0.0 

ERROR2Z=0.0 

IF (N-6) 1000,1000,1030 
gee! 11020 ,1020 ,1010 
LINK=1 


N=I 

IF (X-COEFFT(N)) 1040,1040,1620 
N=I 

GO TO (1130,1060,1190,1230,1080,1090,1180,1330) ,N 
X=Z 

N=2 

TEMP=X 

GO TO 1190 
PRESS=EXP(Z) 
GO TO 1150 

lel =2 

sG=3 
TOLER=.005 

GO TO 1100 
[1=3 

a2 
TOLER=.000005 
H=Y 

LINK=2 

N=1 
X=ALOG(PRESS) 
PLOG=X 

GO TO 1190 
TEMP=Z 

X=TEMP 

=) 
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1160 


1170 
1180 


1190 


1200 
1210 


1220 
1230 
1240 


1250 
1260 


1270 
1280 


1290 
1300 
Hse 
P20 


N=4 

ene) 1p) Ne 

SATLIQ(J)=Z 

JeJ+1 

N=N+1 

he 4J-3) 1190,1190,1170 
IF (LINK) 1640,1640,1300 


3*N-19 


mM) -< >< 


(X-COEFFT(I)) 1200,1200,1210 
Z 

+K 
=X-COEFFT(T) 


NM PH As KRo <M 


nanuvuni unu tt 


1 
I 
G 
0.0 

DO 1220 K=1,10 

I=I+] 

Z=Z*ARG+COEFFT (I) 

GO TO (1140,1070,1050,1160,1160,1160,1410) ,N 

TEMP=Y 

GO TO 1300 

ERROR=H-STEAM( ITI) 

IF (ABS(ERROR)-TOLER) 1640,1250,1250 

IF (ERROR1L) 1260,1280,1260 

ERRORL=ERRORL-ERROR 

IF (ABS(ERROR1)-TOLER) 1290 ,1290,1270 

SLOPE=DELTAL/ERRORL 

DELTAL=SLOPE*ERROR 

ERRORI=ERROR 

TEMP=TEMP+DELTAL 

IF (PRESS-COEFFT(1)) 1320,1320,1310 

IF ((.00001178*PRESS+.09411)*PRESS+382.1-TEMP) 1320,1320,1620 
T= .55555556* TEMP+255. 38223 

TLOG=AL0G(T) 

TAU=1.0/T 

COEFFT (222) =TAU*TAU*186210 .06 

COEFFT(217)=EXP( COEFFT(222)+7.8791476-TLOG) 

COEFFT (221)=162460.*TAU 

COEFFT(220) =126970. *TAU 

COEFFT (219) =EXP (149 .27765-23.*TLOG) 
COEFFT(212)=1.89-COEFFT(217) 

COEFFT (211)=82 .546-COEFFT (221) 
COEFFT(210)=.21828*T-COEFFT(220) 

COEFFT (209)=COEFFT (219 )- .0003635*T 
COEFFT(222)=(COEFFT(222)+COEFFT(222)+1.0) *COEFFT(217) 

COEFFT (217)=COEFFT (212)-COEFFT (222) 

COEFFT(214) =COEFFT(217) /COEFFT(212) 

COEFFT (216) =COEFFT (211) *COEFFT (214)+41 .273-COEFFT (221) 
COEFFT(215) =COEFFT(210) *COEFFT(214) -.5*COEFFT(2 20) 
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1330 


1340 
1350 


COEFFT (214) =COEFFT(209) *COEFFT(214)+1.8461538*COEFFT(219) 
COEFFT (221)=.50*COEFFT (211)-COEFFT (216) 
COEFFT(220)=.25*COEFFT(210)-COEFFT(215) 

COEFFT (219)=.076923077*COEFFT (209 )-COEFFT (214) 

COEFFT(218)=( .0037299965*T+14, 531813) *T+1 7806. 513+472.2493*TLOG 
COEFFT (223)=- .0074599931*T-14 53181 3*TLOG#2.4524439+472 24937 
*4T AU 

P= .068046190*PRESS 

PLOG=ALOG(PRESS) 


SUEFEN(Z13)=COEFPT (212) *1AU 
eoepen(z2o) -COEFFI(213)7P 
mUemet i225) =COERF 1) (226) *COEFF TI (226) 
POPE e 24) COEFF I(225) *COEFF I(225) 
COEFFT (224)=COEFFT (224) *COEFFT (224) *COEFFT (226) 
ome mle 2o) COEFF I(226)*COEFFT(213) 
GOEFFI(213)=P*TAU 
STEAM(1 )=( ( (COEFFT(209) *COEFFT(224)+COEFFT(210) )*COEFFT(225)+ 
SeOerriicl)) )*COEFF I (226)+COEFFT (212 )+4.55504/COEFFT (213) ) 
ae GLo0185 
STEAM(2 )=(( ( (COEFFT(214) *COEFFT(224)+COEFFT(215) )*COEFFT(225)+ 
Sager (2.6) )*COERFT (226 )+COEFFT (217) )*P+COEFFT (218) )*.04355685 
SmeAM(s 1=( 1 ((COEFFI(219) *COEFFT(224)+COEFFI(220) )*COEFFI(225)+ 
eee cree 1) )=COEFFI (220) +COEFF I (222) )*COEFFT (213)+COEFFT (223) 
*+4,55504*P LOG) *(- .02419825) 
QUAL=1 .0 
T=L INK+1 
GO TO (1640,1640,1340,1240,1340,1450) ,I 
X=STEAM(I 1) 
IF (H-X) 1350,1640,1360 
Y=SATLIQ(IT) 
QUAL=(H-Y )/(X-Y) 
X=1.0-QUAL 
STEAM(1 ) =QUAL*STEAM(1 )+X*SATLIQ(1 ) 
STEAM (Jd ) =QUAL*STEAM (Jd )+X*SATLIQ (Jd ) 
STEAM( 11) =H 
IF (LINK-2) 1640,1640 ,1460 


1360 I=(II-1)*11+154 


1370 
1380 


J=1 
X=P LOG 
DO 1380 K=209 ,219 
COEFFT(K)=COEFFT(I) 
Sisal 
X=X-COEFFT(209) 
Y=Y-COEFFT (210) 
vere (209)=(COERFI(211)*X+COEFFT(212) )*X+COEFFT (213) 
COEFFT (210)=(COEFFT (214) *X+COEFFT (215) ) *X+COEFFT (216) 
mOeEr m2 11 )=( (COEFFT(217)*X+COEFFT(218) )*X+COEFFT(219) )*Y 
BSCOEGr b(211 )+COEFFI (219) 
SLOPE=COEFFT(211)+Z 
B21 +COEFF (209) 
GO TO (1390,1430,1530,1540) ,J 
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1390 
1400 
1410 


1420 
1430 


1440 


1450 
1460 
1470 
1480 
1490 
1500 


1510 
fo20 


GN, 


1540 
kote) @ 
1560 


lov’ 
1580 


LINK=3 

IF (Z-TEMP) 1300,1300,1400 
Wee sr 

GO TO 1300 

J=2 

l=22/ 

X=H 

Meeiz-H\) 1420,1370,1370 
1=187 

GO TO 1370 

PED S=oLOPE 

PRESS=EXP(Z) 

PLOG=Z 

IF (1-227) 1520,1520,1440 
LINK=4 

Tl=2 

JJU=3 

GO TO 1120 

DH=H-STEAM (2) 

STEAM(3 )=STEAM(3 )+DH/(TEMP+459.69) 
DS=S-STEAM(3) 
DELTAL=STEAM(3)-ERRORI1 

IF (ERROR1) 1470,1490,1470 
IF (ABS(DELTAL)-.000005) 1500,1500,1480 
DPDS=(PLOG-PLOG1) /DELTAL 
PLOG1=P LOG 

ERRORL=STEAM(3) 

Bee oo-PRESS*(1.0+DS“DPDS) 
IF (LINK-4) 1510,1510,1550 
IF (ABS(DS)-.000005) 1640,1120,1120 
LINK=5 

X=H 

Y=$ 

[=198 

J=3 

fon) 0 1370 

ep = Z 

DTDS=SLOPE 

A=PLOG 

Y=H 

T=165 

J=4 

GO TO 1370 

DTDH=SLOPE 

GO TO 1300 
TEMP=TEMP+DS*DTDS 

IF (ABS(DH)-.005) 1610,1560,1560 
DELTA2=DH-ERROR2 
MeeRRORZ) 157015901570 
IF (ABS(DELTA2)-.005) 1600,1600,1580 
DTDH=(TEMP1-TEMP) /DELTA2 
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£590 
1600 
1610 
1620 
1630 


1640 


TEMP 1=TEMP 

ERROR2=DH 

TEMP=TEMP+DH*DTDK 

a0 1) se 

IF (ABS(DS)-.000005) 1640 ,1600 ,1600 
WRITE (6,1630) NI,XI,YI 

FORMAT (13,2E14.8,13H OUT OF RANGE) 
STOP 

RETURN 

END 
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BRAYTON CYCLE: PRESSURE RATIO=2.00 
WORKING PRESSURE=135.0 

TOTAL PRESSURE DROP=.10 

MASS FLOW (Lbm/s)=444. 37 


STATE TEMP(deg R) SUMMARY:THERMAL EFF=0. 370 
1 540.0 HEATER EFF=0.933 (PHI=0.90) 
2 675.7 SFC(Lbm/HP-hr )=0.572 
3 1668.8 VO TOMA eu £t)=510. 5 
4 2060.0 POWER(HP)=25000.0 
§ 1804.2 WORK FRAC BY STM CYCLE=0.045 
6 Sileieer2 
PP 784.0 
7 Cee 
COMPONENT DPP DHD  VOL(cu ft) 
REGEN £078 =. 744 15.5 EFF=0.88 
TURBINE 10.0 EFF(P)=0.91 
COOLER 2007 B20 139.8 
COMPRESSOR 29.3 EFF(P)=0.88 
HEATER 0005 
BOILER O14 EFF=0.94 
HEATER CYCLE: PRESSURE RATIO = 4.79 
TOTAL PRESSURE DROP=.43 
AIR FLOW(Lbm/s) =47.84 
FUEL FLOW(Lbm/s)= 3.97 
STATE TEMP(deg R) 
th 535.0 
2h 315)5) 52 
3h 77 a 
Wh 5054.5 
5h 1768.8 
6h 945.0 
7h 760.0 
COMPONENT DPP DHD  vVOL(cu ft) 
RECUP  .321 650 a : EFF=0.96 
TURBINE 
HEATER S02 > 310 Pes ‘i 
COMPRESSOR 16.0 EFF(P)=0.85 


RANKINE CYCLE 


STEAM FLOW(Lbm/s)= 
.0067 


FLOW RATIO=0 


SATURATION PRESSURE= 
2.96 


Oma 


STATE TEMP( deg R) PRESSURE( psi) 
1s 540.0 Or 5 
2S S566 SO 
SF 749.0 56.7 
3s TRS fe Oo" 
Us 540.0 0.5 QUAL=0.844 
COMPONENT base che iat) 
BOILER 162.3 
TURBINE in EFF(IS )=0.88 
CONDENSER 68 .3 
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BRAYTON CYCLE: PRESSURE RATIO= 3.00 
WORKING PRESSURE =135.0 
TOTAL PRESSURE DROP=.10 
MASS FLOW(Lbm/s)=287.63 


STATE TEMP(deg R) SUMMARY: THERMAL EFF=0.432 
i 540.0 HEATER EFF=0.933 (PHI=0.90) 
2 769.8 SFC(Lbm/HP-hr) =0.490 
3 LS) S16) & fe ViGion LO eC leert) = 586. 1 
4 Z0e0.0 POWER (HP) =25000.0 
5 1643.6 WORK FRAC BY STM CYCLE=0.044 
6 874.7 
PP 841.2 
7 eats 
COMPONENT DPP DHD  VOL(cu ft) 
REGEN .050 .837 11.5 EFF=0.88 
TURBINE 5.3 EFF(P)=0.91 
COOLER .024 .352 88.9 
COMPRESSOR 20.9 EFF(P)=0.88 
HEATER 0004 
BOILER .026 EFF=0.94 


HEATER CYCLE: PRESSUPE RATIO= 4.79 
TOTAL PRESSURE DROP=.43 

AIR FLOW(Lbm/s)= 40.99 

FUEL FLOW(Lbm/s)= 3.40 


STATE TEMP(deg R) 
ih SiS ne, 
2h 835.2 
3h 1586.5 
Lh 4932.9 
Sh 1638.8 
6h 945.0 
OG. 70.0 
COMPONENT DPP DHD- VOL(cu ft) 

RECUP 10 650 2.2 EFPF=0.93 
TURBINE 0.9 EFF(P)=0.89 
HEATER 6022 .358 ee) Ble. 

COMPRESSOR P2ny7) SRE CE )A0.65 


RANKINE CYCLE: SATURATION PRESSURE= 128.0 
STEAM FLOW(Lbm/s)= 2.48 
FLOW RATIO=0.0086 


SUD TEMP(deg R) PRESSURE (psi) 
1s 540.0 Oss 
2S 555.0 eZ 0 
Sr 806.2 ZG O 
3s 839.7 125.0 
Ls 540.0 0.5 QUAL=0.819 
COMPONENT VOL(cu ft) 
BOILER 638.3 
TURBINE 3.4 EFF(IS)=0.88 
CONDENSER Som) 
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BRAYTON CYCLE: 


WORKING PRESSURE = 1 
TOTAL PRESSURE DROP=.10 
MASS FLOW (Lbm/s)=239.91 


STATE 


OWN mu) Ne 


‘aaa 


7. 
COMPONENT 
REGEN 
TURBINE 
COOLER 
COMPRESSOR 
HEATER 
BOILER 


BEATER CYCLE: 


PRESSURE RATIO= 4.00 


35.0 


TOTAL PRESSURE DROP=.43 
AIR FLOW(Lbm/s)=39.59 
FUEL FLOW(Lbm/s)= 3.29 


STATE 
Lh 
ch 
3n 
4h 
5h 
6h 
7n 

COMPONENT 

RECUP 

TURBINE 
HEATER 
COMPRESSOR 


RANKINE CYCLE: 


TEMP( deg R) SUMMARY: THERMAL EFF=0.447 
540.0 HEATER EFF=0.933 (PHI=0.90) 
843.9 SFC(Lbm/HP-hr )=0.473 

1453.8 VOMeLORAnVcur mt) =) 352.2 
2060.0 POWER(HP)=25000.0 
1536.9 WORK FRAC BY STM CYCLE=0.063 
O27. 0 

876.8 

860.6 
DPP DHD VOoL(cu ft) 

FOZ 2950 13.2 EFF=0.88 

3.6 EFF(P)=0.91 
6032 3.444 89.8 
18.4 EFF(P)=0.88 

.0003 

O46 EFF=0.94 
PRESSURE RATIO= 4.79 

TEMP(deg R) 

B55 50 

835.2 
1494.5 
4855.0 

SO SIGs 

945.0 

760. 0 
DPP DHD VOibeeu, 2) 

.299 .650 2.1 EFF=0.92 

0.8 EFF(P)=0.89 
.033 .406 dace, 


STEAM FLOW(Lbm/s) = 


FLOW RATIO=O. 


STATE 
1s 
2S 
SF 
3s 
Is 
COMPONENT 
BOILER 
TURBINE 
CONDENSER 


Ope S15) 


12.0 EFF(P)=0.85 


SATURATION PRESSURE= 199.9 


3.19 


some ee R) PRESSURE( psi) 
Om 


540.0 
oe 0 199.9 
841.8 199.9 
892.0 199.9 
540.0 O25.) GUAT =0 .510 
VOlCeu tt ) 
119.0 
Speer (iS }=0.08 


70.5 
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BRAYTON CYCLE: PRESSURE RATIO= 5.00 
WORKING PRESSURE=135.0 

TOTAL PRESSURE DROP=.10 

MASS FLOW (Lbm/s)=212. 34 


STATE TEMP(deg R) SUMMARY: THERMAL EFF=0.460 
1 540.0 HEATER EFF=0.933 (PHI=0.90) 
2 905.8 SFC(Lbm/HP-hr)=0.460 
2 1391.8 VOL TOTAL(cu ft)=596.1 
2060.0 POWER( HP) =25000.0 
5 1458.1 WORK FRAC BY STM CYCLE=0.111 
6 See 
PP 876.8 
e 846.9 
COMPONENT DPP DHD  VOL(cu ft) 
REGEN Oli, Cane S0 15.7 EFF=0.88 
TURBINE 2.7 EFF(P)=0.91 
COOLER 022 3530 LOvac 
COMPRESSOR 16.7 EFF(P)=0.88 
HEATER 0003 
BOILER 068 EFF=0.94 


HEATER CYCLE: PRESSURE RATIO= 4.79 
TOTAL PRESSURE DROP=.43 

AIR FLOW(Lbm/s)= 38.49 

FUEL FLOW(Lbm/s)= 3.19 


STATE TEMP(deg R) 


th SS ene, 
2n 835.2 
3h 1427.4 
Ln 4798.9 
Sh 1491.8 
6h 945.0 
7h PASO) 10 
COMPONENT DPP  ODHD VOL(cu ft) 

RECUP .289 .650 2 EFF=0.90 
TURBINE 0.8 EFF(P)=0.89 
HEATER O44 (454 12.4 

COMPRESSOR 11.5 EFF(P)=0.85 


RANKINE CYCLE: SATURATION PRESSURE= 200.0 
STEAM FLOW(Lbm/s)= 5.21 
FLOW RATIO=0.0245 


STATE TEMP(deg R) PRESS URE( psi) 
On5 


is 540.0 
2s 555.60 200.0 
Sieh 841.8 ZOO 0 
3S 937.1 200,00 
lbs 540.0 0.5 QUAL=0.826 
COMPONENT voL(cu ft) 
BOILER 298.1 
TURBINE 10.9 EFF(IS)=0.88 
CONDENSER 117.4 
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BRAYTON CYCLE: PRESSURE RATIO = 6.00 
WORKING PRESSURE=135.0 
TOTAL PRESSURE DROP=.10 
MASS FLOW(Lbm/s)=197.00 


STATE TEMP( deg R) SUMMARY: THERMAL EFF=0.464 


i] 540.0 HEATER EFF=0.933 (PH1I=0.90) 
2 959.4 SFC(Lbm/HP-hr )=0.456 
3 1343.8 VOL TOTAL(cu ft) =879.0 
h 2060.0 POWER(HP)=25000.0 
5 1396.2 WORK FRAC BY STM CYCLE=0.151 
6 1011.8 
PP 876.8 
SSS, 
COMPONENT DPP DHD  VOL(cu ft) 
REGEN NO 6 o50 14.4 EFF=0.88 
TURBINE 2.2 EFF(P)=0.91 
COOLER LOWS. Wnts AKC) 9 
COMPRESSOR 15.8 EFF(P)=0.88 
HEATER .0003 
BOILER 074 EFF=0 .94 


HEATER CYCLE: PRESSURE RATIO= 4.79 
TOTAL PRESSURE DROP=.43 

AIR FLOW(Lbm/s)=38.17 

FUEL FLOW(Lbm/s)= 3.17 


STATE TEMP(deg R) 


th 535.0 
2h 835.2 
3h LAG ie 
Lh 4755.8 
Sh 1443.8 
6h 945.0 
7h 760.0 
COMPONENT DEP SHB! VOL(cu ft) 

RECUP 2278) 2050 al EFF=0.89 
TURBINE 0.8 EFF(P)=0.89 
HEATER .055 .502 9.5 

COMPRESSOR 11.4 EFF(P)=0.85 


RANKINE CYCLE: SATURATION PRESSURE= 200.0 
STEAM FLOW(Lbm/s)= 6.70 
FLOW RATIO=0.0340 


SEND TEMP(deg R) PRESSURE( psi) 


is 540.0 0.5 
2s Do) 20 200.0 
Syl 841.8 200.0 
3s 976.8 200.0 
Ls 540.0 O25 GUAL-0.6 36 
COMPONENT VOmCeu ft) 
BOILER 523.6 
TURBINE 16.2 EFF(IS)=0.88 
CONDENSER ILS 6 
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Beer ON CYCLE: PRESSURE RATIO=7.00 
WORKING PRESSURE=135.0 
metas PRESSURE DROP=.10 
MASS Flow(Lbm/s) =187.48 


STATE TEMP(deg R) SUMMARY: THERMAL EFF=0.463 


1 540.0 HEATER EFF=0.933 (PHI=0.90) 
2 100720 SFC(Lbm/HP=hr) =0.457 
3 1304.9 VOL TOTAL(cu ft) =1195.7 
4 2060.0 POWER(HP) =25000.0 
5 1345.5 WORK FRAC BY STM CYCLE=0.189 
6 1047.6 
PP 876.8 
7 825.2 
COMPONENT DPP DHD VOL(cu ft) 
REGEN sOmOr =, 9150 13.6 EFF=0.88 
TURBINE 1.9 EFF(P)=0.91 
COOLER "Oles 2720 ILS AS, 
COMPRESSOR jag DE Cr )i=0.36 
HEATER . 0004 
BOILER .078 EFF=0.94 
HEATER CYCLE: PRESSURE RATIO= 4.79 
TOTAL PRESSURE DROP=.43 
AIR FLOW(Lbm/s)= 38.22 
FUEL FLOW(Lbm/s)= 3.17 
STATE TEMP(deg R) 
th 535.0 
2h 835.2 
3n ooo 
Lh LPAI ae Pe 
Sh 1404.9 
6h 945.0 
7h 760.0 
COMPONENT DPP DHD- VOL(cu ft) 
RECUP 720" 2.050 2.1 EFF=0.87 
TURBINE 0.8 EFF(P)=0.89 
HEATER .066 .550 fA) 
COMPRESSOR 11.4 EFF(P)=0.85 


RANKINE CYCLE: SATURATION PRESSURE= 200.0 
STEAM FLOW(Lbm/s)= 7.92 
FLOW RATIO=0.0422 


STATE TEMP(deg R) PRESS URE( psi) 


is 540.0 0.5 
2s 555 «0 200.0 
Sue 841.8 ZUG. 0 
3s LCA ate 208.0 
Ls 540.0 0.5 QUAL=0.849 
COMPONENT Vomcu £%) 
BOILER 782.4 
TURBINE Ae ohhh 15 )=0.88 
CONDENSER 183.6 
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BRAYTON CYCLE: PRESSURE RATIO= 8.00 
WORKING PRESSURE= 135.0 


TOTAL PRESSURE DROP= .10 
MASS FLOW(Lbm/s) = 181.22 
STATE TEMP(deg R) SUMMARY 
1 540.0 
2 1049.8 
3 1272e5 
Ly 2060.0 
5 1302.9 
6 1080.2 
PP 876.8 
7 SiGe 3 
COMPONENT DPP  ODHD VoOmGeu £% 
REGEN “O10 O50 sles yeel 
TURBINE 6 
COOLER FO) OS Wet, fal 
COMPRESSOR 15.1 
HEATER .0004 
BOILER 080 


HEATER CYCLE: PRESSURE RATIO= 4.79 
TOTAL PRESSURE DROP=.43 

AIR FLOW(Lbm/s) =38.45 

FUEL FLOW(Lbm/s)= 3.19 


:THERMAL EFF=0.460 

HEATER EFF=0.933 (PHI=0.90) 
SFC(Lbm/HP-hr) =0.459 

VOL TOTAL(cu ft)=1550.7 
POWER (HP)=25000.0 

WORK FRAC BY STM CYCLE=0.225 


) 

EFF=0.88 
EFF(P)=0.91 
EFF(P)=0.88 


EFF=0.94 


STATE  TEMP(deg R) 
Lh 535.0 
2h 835.2 
3h 1298.2 
Lh 4692.5 
5h 1a7e.> 
6h 945.0 
7h 760.0 
COMPONENT DPP DHD VOL(cu ft) 

RECUP .256 .650 Zen Err. 66 
TURBINE 0.8 EFF(P)=0.89 
HEATER 076 .598 6.5 

COMPRESSOR 11.5 EFF(P)=0.85 
RANKINE CYCLE: SATURATION PRESSURE=200.0 
STEAM FLOW(Lbm/s)= 8.97 
FLOW RATIO=0.0495 

SAIS: TEMP(deg R) PRESS URE( psi) 

1s 540.0 Om5 
2s bins) 18 ZO0n0 
Sug 841.8 200.0 
es 1045.2 200.0 
Ss 540.0 0.5 QUAL=0.859 
COMPONENT VO leceumsct ) 
BOILER GS) 
TURBINE Zoomer (is }=0.68 
CONDENSER Z21LOee 
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BRAYTON CYCLE: 


STATE 


Ig 


“NJ'IG OWN Fly NP 


COMPONENT 
REGEN 
TURBINE 
COOLER 
COMPRESSOR 
HEATER 
BOILER 


HEATER CYCLE: 


PRESSURE RATIO= 9.00 
WORKING PRESSURE=135.0 

TOTAL PRESSURE DROP=.10 

MASS FLOW(Lbm/s)=176.95 


TEMP(deg R) 


540.0 
1089.0 
1245.0 
2060.0 
2OCn2 
WIL LO) 62 

876.8 

SOowc 
DPP 
2010 


.0O10 


0004 
080 


DHD 


SUMMARY: THERMAL EFF=0.456 
HEATER EFF=0.933 (PHI=0.90) 
SFC(Lbm/HP-hr)=0.464 
VOL TOTAL(cu ft)=1950.3 
POWER (HP )=25000.0 
WORK FRAC BY STM CYCLE=0.260 


VOmcenm ft ) 
L217 Feh=0.66 

ime PPPCP)=0.901 

1) 

15.0 


-950 


904 
EFF(P)=0.88 


EFF=0.94 


PRESSURE RATIO= 4.79 


TOTAL PRESSURE DROP = 
AIR FLOW(Lbm/s)= 38.80 
FUEL FLOW(Lbm/s)= 3.22 


STATE 
ih 
Za 
3h 
Lh 


COMPONENT 
ree UP 
TUR BINE 
HEATER 
COMPRESSOR 


RANKINE 


STATE 
1s 
2S 
oF 
3s 
ls 
COMPONENT 
BOILER 
TURBINE 
CONDENSER 


C1 CLE: 


DEY 
245 


087 


DHD 


43 


TEMP(deg R) 


VORGcu ft) 
Bae ee =05 65 
0.8 EFF(P)=0.89 


.650 


646 Or 
11.7 EFF(P)=0.85 


SATURATION PRESSURE= 200.00 
STEAM FLOW(Lbm/s)= 9.90 
FLOW RATIO=0.0560 


TEMP( deg R) PRESSURE( psi) 

540.0 : 

So. 0 2000 

841.8 ZO 19 
Me 7 5-2 ZOOmO 

540.0 O75 GURE=0.867 
VOL(cu ft) 

1439.0 

Zon seer (is) =0.08 
24.5 
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BRAYTON CYCLE: PRESSURE RATIO=10.00 
WORKING PRESSURE=1 35.0 

TOTAL PRESSURE DROP=.10 

MASS FLOW (Lbm/s)=174.00 


STATE TEMP(deg R) SUMMARY: THERMAL EFF=0.451 


1 540.0 HEATER EFF=0.933 (PHI=0.90) 
2 1125.1 SFC(Lbm/HP-hr)=0.469 
? LZ 2sL sal VOL TOTAL(cu ft) =2403.0 
2060.0 POWER (HP) =25000.0 
5 1234.2 WORK FRAC BY STM CYCLE=0.296 
6 1138.1 
PP See) 7 
7 800.9 
COMPONENT DPP DHD VOL(cu ft) 
REGEN satel aisiGre 12.4 EFF=0.88 
TURBINE 1.3 EFF(P)=0.91 
COOLER 2008 .996 232.8 
COMPRESSOR 14.9 EFF(P)=0.88 
HEATER 0004 
BOILER .082 EFF=0.94 


HEATER CYCLE: PRESSURE RATIO= 4.79 
TOTAL PRESSURE DROP=.43 

AIR FLOW(Lbm/s)= 39.22 

FUEL FLOW(Lbm/s)= 3.25 


STATE TEMP(deg R) 


2h 835.2 
3h 1242.5 
4h 4647.3 
Sh ied. 1 
6h 945.0 
7h 760.0 
COMPONENT DPP DHD  VOL(cu ft) 

RECUP 1235 ©6650 2.2 EFF=0.84 
TURBINE 0.8 EFF(P)=0.89 
HEATER .098 .694 5.0 

COMPRESSOR io EFPR(P)=0.05 


RANKINE CYCLE: SATURATION PRESSURE= 199.9 
STEAM FLOW(Lbm/s)= 10.76 
FLOW RATIO=0.0618 


LATE TEMP( deg R) PRESSURE( psi) 
540.0 


1s , 
2S Sas ae 199.9 
SF 841.7 199.9 
35 POS. 1 sRelodne 
is 540.0 0.5 QUAL=0.875 
COMPONENT vOL(cu ft) 
BOILER 1830.6 
TURBINE 34.0 EFF(IS)=0.88 
CONDENSER 257.1 


88 





Thesis 207056 
| BB75 Burgamy 
e@2 Preliminary design of 
an alternative fuels 
combined cycle propul- 
sion plant for naval 
ship applications. 


Thesis 207956 
B875 Burgamy 
Che Preliminary design o° 
an alternative fuels 
combined cycle propul- 
sion plant for naval 
ship applications. 








thesB875 
Preliminary design of an alternative fue 


HAI 


3 2768 002 07933 7 
DUDLEY KNOX LIBRARY 





